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ON EXPLICIT SOLUTIONS TO THE STATIONARY AXISYMMETRIC 
EINSTEIN-MAXWELL EQUATIONS DESCRIBING DUST DISKS 

C. KLEIN 


Abstract. We review explicit solutions to the stationary axisymmetric Einstein-Maxwell 
equations which can be interpreted as disks of charged dust. The disks of finite or infi¬ 
nite extension are infinitesimally thin and constitute a surface layer at the boundary of an 
electro-vacuum. The Einstein-Maxwell equations in the presence of one Killing vector are 
obtained by using a projection formalism. This leads to equations for three-dimensional 
gravity where the matter is given by a SU( 2, 1)/S[U(1, 1) x U( 1)] nonlinear sigma 
model. The SU (2,1) invariance of the stationary Einstein-Maxwell equations can be used 
to construct solutions for the electro-vacuum from solutions to the pure vacuum case via 
a so-called Harrison transformation. It is shown that the corresponding solutions will al¬ 
ways have a non-vanishing total charge and a gyromagnetic ratio of 2. Since the vacuum 
and the electro-vacuum equations in the stationary axisymmetric case are completely in¬ 
tegrate, large classes of solutions can be constructed with techniques from the theory of 
solitons. The richest class of physically interesting solutions to the pure vacuum case due 
to Korotkin is given in terms of hyperelliptic theta functions. The Harrison transformed 
hyperelliptic solutions are discussed. As a concrete example we study the transformation 
of a family of counter-rotating dust disks. To obtain algebro-geometric solutions with van¬ 
ishing total charge which are of astrophysical relevance, three-sheeted surfaces have to be 
considered. 

The matter in the disk is discussed following Bicak et al. We review the ‘cut and glue’ 
technique where a strip is removed from an explicitly known spacetime and where the re¬ 
mainder is glued together after displacement. The discontinuities of the normal derivatives 
of the metric at the glueing hypersurface lead to infinite disks. If the energy conditions 
are satisfied and if the pressure is positive, the disks can be interpreted in the vacuum case 
as made up of two components of counter-rotating dust moving on geodesics. In electro¬ 
vacuum the condition of geodesic movement is replaced by electro-geodesic movement. 
As an example we discuss a class of Harrison-transformed hyperelliptic solutions. The 
range of parameters is identified where an interpretation of the matter in the disk in terms 
of electro-dust can be given. 


1. Introduction 

Electromagnetic fields, especially magnetic fields play a role in astrophysics in the con¬ 
text of neutron stars, white dwarfs and galaxy formation. A complete relativistic under¬ 
standing of such situations requires studying the coupled Einstein-Maxwell equations. Of 
special interest are stationary axisymmetric situations since isolated matter configurations 
in thermodynamical equilibrium belong within relativity to this class, see (Dll' Since 
the stationary axisymmetric Einstein-Maxwell equations in vacuum in the form of Ernst 
0 are completely integrable (see 0), powerful solution generating techniques from the 
theory of solitons are at hand to obtain physically interesting solutions. But the equations 
in the matter region — which is generally approximated as an ideal fluid — do not seem 
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to be integrable. This makes it difficult to find global solutions which hold both in a three- 
dimensionally extended matter region and in vacuum. 

To obtain global solutions to the Einstein-Maxwell equations with these methods, one 
is thus limited to two-dimensionally extended matter distributions, i.e. surface layers. In¬ 
finitesimally thin disks have been discussed in Newtonian astrophysics as models for cer¬ 
tain galaxies, see 0, and as models for the matter in accretion disks around black holes, 
see 0 and references therein. In this case the equations in the matter region reduce to 
ordinary differential equations the solutions of which determine boundary data for the vac¬ 
uum equations. Alternatively surface layers can be obtained by ‘cut and glue’ techniques. 
Bicak, Lynden-Bell and Katz |7 8? studied static spacetimes, from which they removed 
a strip and glued the remainder together. The non-continuous normal derivatives at the 
gluing plane lead to a (5-type energy-momentum tensor which can be interpreted as an in¬ 
finitely extended disk made up of counter-rotating dust. This method was extended in j9j 
to generate disk sources of the Kerr-metric. With the same techniques, disk sources for 
Kerr-Newman metrics m, static axisymmetric spacetimes with magnetic fields CD and 
conformastationary metrics ca were given. 

Counter-rotating disks are discussed in astrophysics as models for certain SO and Sa 
galaxies (see CD and references given therein). These galaxies show counter-rotating 
matter components and are believed to be the consequence of the merger of galaxies. Re¬ 
cent investigations have shown that there is a large number of galaxies (see GZI, the first 
was NGC 4550 in Virgo) which show counter-rotating streams in the disk with up to 50 % 
counter-rotation. 

By construction all disks due to ‘cut and glue’ techniques have an infinite extension but 
finite mass since the mass of the spacetime is not changed by the method. The matter in 
the disks can be interpreted as a two-dimensional fluid with a purely azimuthal pressure. 
If the energy conditions in the disk are satisfied and if the pressure in the disk is positive, 
the matter can alternatively be interpreted as consisting of two counter-rotating streams of 
pressureless matter, so-called dust. In the pure vacuum case this is best done by introducing 
observers rotating with the disk in a way that the energy-momentum tensor is diagonal for 
them. It can be shown that the corresponding dust streams move on geodesics of the 
inner geometry of the disk. In the electro-vacuum case, the corresponding condition on 
the matter is motion on electro-geodesics, i.e. solutions to the geodesic equation in the 
presence of a Lorentz force. It was shown in ED that this is a more restrictive condition 
than in the pure vacuum case. 

For vacuum so-called Riemann-Hilbert techniques (see 1141 03 s 161 and references 
therein) were used to generate solutions for disks of finite extension. Explicit metrics 
could be given in terms of theta functions on hyperelliptic Riemann surfaces fl4l . Since 
these two-sheeted surfaces are a generalization of the well-known elliptic surfaces, a pow¬ 
erful theory is at hand to treat hyperelliptic functions. The main advantages of this class 
are that it is very rich (‘solitonic’ solutions as the Kerr solution for a rotating black hole 
are contained as limiting cases), and that a whole subclass with physically interesting so¬ 
lutions could be identified in 1 1711181 by studying the analyticity properties of the solution. 
The corresponding solutions for the electro-vacuum are given on three-sheeted surfaces 
O which are mathematically less well understood. Recent progress in this context was 
made by Korotkin EO by considering the Riemann-Hilbert problem on multi-sheeted cov¬ 
erings of the complex plane. The solutions to this problem which can be used to solve the 
Einstein-Maxwell equations are again given in terms of theta functions. 
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Until now there are, however, no explicit examples for physically realistic disk solu¬ 
tions on three-sheeted Riemann surfaces. Therefore an intermediate step was taken in 1201 
where hyperelliptic solutions with charge were studied. These solutions were obtained by 
exploiting the SU( 2,1) invariance of the stationary Einstein-Maxwell equations (see 1211 
t0 1251 1. Using a so-called Harrison transformation 1261 . one can generate solutions with 
charge from solutions to the pure vacuum equations. A remarkable property of the cor¬ 
responding spacetimes is the fact that their gyromagnetic ratio is always identical to 2 as 
in the case of the Kerr-Newman black holes. This is of interest in the context of claims 
in t27l f28l that this property is a hint on a deep connection between general relativity 
and relativistic quantum mechanics. By studying the asymptotics of Harrison-transformed 
pure vacuum solutions, it was shown that the thus obtained solutions will always have a 
non-vanishing total charge which limits their astrophysical relevance since charges seem 
to neutralize in our universe. This is a hint that astrophysically interesting solutions with¬ 
out total charge, but non-vanishing magnetic fields in terms of theta functions are to be 
expected only on three-sheeted surfaces. 

This paper is is organized as follows: In section 2 the Newtonian case is studied for 
illustration. The ‘cut and glue’ techniques of Bicak et al. 01 and disks of finite extension 
are presented. In section 3 we consider the Einstein-Maxwell equations in the presence 
of a single Killing vector. Using a projection formalism | 29i 301 , we perform a standard 
dimensional reduction of the equations. It is shown that the stationary Einstein-Maxwell 
equations are equivalent to three-dimensional gravity with a SU (2,1 )/S[U (1,1) x (7(1)] 
sigma model as matter. Using complex notation, one can introduce Ernst potentials 0. 
We study a gauge invariant formulation of the SU(2, 1) matrix of the sigma model and 
the related transformations of the solutions. The Harrison transformation is presented and 
discussed for simple examples. The asymptotic behavior of asymptotically flat solutions 
is studied. The stationary axisymmetric case is shown to be completely integrable. In 
section 4 we recall basic facts on the stationary axisymmetric pure vacuum case and on 
hyperelliptic disk solutions. Using the Harrison transformation on a family of counter¬ 
rotating disk solutions [3J, 32]|, we obtain the complete transformed metric and discuss 
interesting limiting cases. The discussion of the energy-momentum tensor using Israel’s 
junction conditions (33| is presented in section 5. The case of the Harrison transformed 
counter-rotating disk is studied as an example. The range of the physical parameters is 
given where the matter in the disk can be interpreted as electro-dust. In section 6 we 
summarize recent results by Korotkin ED on solutions to the Riemann-Hilbert problem 
on multi-sheeted Riemann surfaces in terms of Szego kernels and solutions to the Einstein- 
Maxwell equations. In section 7 we add some concluding remarks. 

2. Newtonian dust disks 

To illustrate the basic concepts used in the following sections, we will briefly recall 
some facts on Newtonian dust disks. In Newtonian theory, gravitation is described by a 
scalar potential U which is a solution to the Laplace equation in the vacuum region. The 
units in this article are chosen in a way that the Newtonian gravitational constant, the 
dielectric constant and the velocity of light are equal to 1. We use cylindrical coordinates 
p, ( and cj) and place the disk made up of a pressureless two-dimensional ideal fluid with 
radius po in the equatorial plane ( = 0. In Newtonian theory stationary perfect fluid 
solutions and thus also the here considered disks are known to be equatorially symmetric. 

Since we concentrate on dust disks, i.e. pressureless matter, the only force to compen¬ 
sate gravitational attraction in the disk is the centrifugal force. This leads in the disk to 
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(here and in the following f, x = 

(i) u, p = n 2 (p) p , 

where fl(p) is the angular velocity of the dust at radius p. Since all terms in 0 are qua¬ 
dratic in there are no effects due to the sign of the angular velocity. The absence of these 
so-called gravitomagnetic effects in Newtonian theory implies that disks with counter¬ 
rotating components will behave with respect to gravity exactly as disks made up of only 
one component. We will therefore only consider the case of one component in this section. 
Integrating o we get the boundary data U(p,0) with an integration constant Uq = U(0, 0) 
which is related to the central redshift in the relativistic case. 

To find the Newtonian solution for a given rotation law f2(p), we have to construct 
a solution to the Laplace equation which is everywhere regular except at the disk where 
it has to take the boundary data 0 . At the disk the normal derivatives of the potential 
will have a jump since the disk is a surface layer. Notice that one only has to solve the 
vacuum equations since the two-dimensional matter distribution merely leads to boundary 
conditions for the Laplace equation. In the Newtonian setting one thus has to determine 
the density for a given rotation law or vice versa, a well known problem (see e.g. 0 and 
references therein) for Newtonian dust disks. 

There are several ways to construct Newtonian dust disks. We will only outline two 
possibilities which can be used with some modifications also in the relativistic case. 


2.1. ‘Cut and glue’-techniques. One way to construct Newtonian disks is to start with a 
known equatorially symmetric solution to the Laplace equation, for instance the solution 
for a point mass m, 


( 2 ) 


U=- 


m 

vV + C 2 ' 


Then a strip of width 2£o is cut out of the space symmetrically to the equatorial plane. The 
solutions for positive and negative £ are displace in ^-direction by ±£o and glued together 
at the equatorial plane. The discontinuity at this plane leads to a surface layer of infinite 
extension. 

In the example of a point mass considered by Kuzmin d and Toomre G3, this leads 
to the solution 


(3) 


vV + aci+Co) 2 ’ 


The surface density ad at the equatorial plane is just given by 2nad = (0 + ), in the 

example 


(4) 


u c ( o + ) 


m( 0 

{P 2 + Q)i~ 


By construction the spacetime has finite mass m since the asymptotics have not been 
changed by the procedure of cutting and glueing. The disk is infinite, but the mass density 
decreases because of @ as p 3 . Since the mass density is positive, the matter in the disk 
can be interpreted as a disk of dust. The angular velocity as defined in 0 is in this example 
given by 


n 2 = 


in 


O 2 + CoP 


(5) 
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Asymptotically the angular velocity satisfies the Kepler relation f l 2 p 3 = m for test parti¬ 
cles. Thus the matter in the disk behaves for large distances from the center as free par¬ 
ticles. Though the gravitational field □ can be seen as generated by the self-gravitating 
matter in the disk, particles in large distances from the center can be viewed as test particles 
since the the density tends to zero. The so constructed disks thus have physically accept¬ 
able properties: a positive mass density and a finite mass. Though they have an infinite 
extension, the density decreases rapidly for p —> oo. 

Since the Laplace equation is linear, arbitrary linear combinations of potentials of the 
form 0 will always leads to infinite disks. Evans and de Zeeuw f36| considered disk 
potentials of the form 


( 6 ) 


TT= [ v< y e ) de 

J VP 2 + ( ICI+e ) 2 


which leads to the general classical disk formula. Bicak, Lynden-Bell and Katz Q used 
this technique to generate disk solutions to the static axisymmetric Einstein equations. 


2.2. Disks of finite extension. To generate disks of finite extension, we use an approach 
which can be generalized to some extent to the relativistic case. The resulting expression 
will be shown to be equivalent to the Poisson integral for a distributional density. We put 
p 0 = 1 without loss of generality (we are only considering disks of finite non-zero radius) 
and obtain U as the solution of a so-called Riemann-Hilbert problem (see e.g. E6) and 
references given therein). The solution can be written in the form 


(7) 


U(p,Q 


1 f In G{K)dK 

J t ^K - C ) 2 + p 2 ’ 


where In G £ C' 1, “(T) (Holder continuous on T) and where T is the covering of the imag¬ 
inary axis in the upper sheet of C between —i and i; C is the Riemann surface of genus 0 
given by the algebraic relation Pq(K) = ( K — £) 2 + p 2 . The function G has to be subject 
to the conditions G(K) = G(K) and G(—K) = G(K). 

It may be checked by direct calculation that U in 0 is a solution to the Laplace equa¬ 
tion except at the disk. The reality condition on G leads to a real potential, whereas the 
symmetry condition with respect to the involution I\ —> -K leads to equatorial symmetry. 
The occurrence of the logarithm in 0 is due to the Riemann-Hilbert problem with the help 
of which the solution to the Laplace equation was constructed, see e.g. M- 

The function In G is determined by the boundary data U(p, 0) or the energy density a,i 
of the dust via 


( 8 ) 


In G(t) 


4 \U 0 +1 


U P (p)dp \ 

V^7J 


or 

,9, ,» G(t ) = 4/‘A rfp 

respectively where t = —iK. This can be seen in the following way: 

At the disk the potential takes due to the equatorial symmetry the boundary values 


( 10 ) 


u(p, 0) 


irm# 

Jo \Jp 2 - f 2 
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and 


( 11 ) 


v ^—hS. 


^(lnGCf)) 


dt. 


V t2 ~ P 2 

Both equations constitute integral equations for the ‘jump data’ InG of the Riemann- 
Hilbert problem if the respective left-hand side is known. The equations ( II 01 1 and flD 
are both Abelian integral equations and can be solved in terms of quadratures, i.e. <0 and 

0 . 

To show the regularity of the potential U, we prove that the integral 0 is identical to 
the Poisson integral for a distributional density which reads at the disk 
( 12 ) 


U(p) = -2 [ a d {p')p'dp' [ 
Jo Jo 


d(j) 


0 \j{p + p') 2 - 4/9(0' COS (j) 


= -4 


a d (p')p'dp‘ 


, , K(k(p , p')) 
P + P' ’ 


where k{p,p') = 2y / pp 7 /(p + p') and where K(k) is the complete elliptic integral of 
the first kind. Eliminating InG in (1101 via 0 we obtain after interchange of the order of 
integration 


03) 




which is identical to ( fill since K(2y/k/(l + k)) = (1 + k)K(k). Thus the integral 
0 has the properties known from the Poisson integral: it is a solution to the Laplace 
equation which is everywhere regular except at the disk where the normal derivatives are 
discontinuous. 

We note that it is possible in the Newtonian case to solve the boundary value problem 
purely locally at the disk. The regularity properties of the Poisson integral then ensure 
global regularity of the solution except at the disk. Such a purely local treatment will not 
be possible in the relativistic case. 

The above considerations make clear that one cannot prescribe both U at the disk (and 
thus the rotation law) and the density independently. This just reflects the fact that the 
Laplace equation is an elliptic equation for which Cauchy problems are ill-posed. If In G 
is determined by either © or 0 for given rotation law or density, expression 0 gives the 
analytic continuation of the boundary data to the whole spacetime. In case we prescribe 
the angular velocity, the constant Uq is determined by the condition lnG(i) = 0 which 
excludes a ring singularity at the rim of the disk. Lor rigid rotation (f l = const), we get 
e.g. 

(14) In G(t) = 4fl 2 (r 2 + 1), 


which leads with 0 to the well-known Maclaurin disk. 


2.3. Charged static dust disks. In a Newtonian theory gravity and electromagnetism de¬ 
couple. Disks of pressureless charged matter will lead to a gravitational potential U as 
above and to electric and magnetic fields. Lor disks with only one component of dust, 
there will be necessarily a magnetic field due to the rotating charges. Static solutions are 
possible if two streams of charged dust with equal densities are exactly counter-rotating. In 
this case the magnetic field vanishes since the effects of both streams concerning magnetic 
fields just compensate. This trick was used by Morgan and Morgan ED within general rel¬ 
ativity to describe static disks: in the case of two identical streams of counter-rotating parti¬ 
cles the so-called ‘gravitomagnetic’ effects of relativity cancel, and the resulting spacetime 
is static. 







CHARGED DUST DISKS 


7 


In the electrostatic case, the electric potential is also a solution to the Laplace equation. 
If we assume the electrical density a e to be proportional to the matter density ad in the case 
of a dust disk (both densities are surface densities), i.e. a e = Qad , we get from Newton’s 
law Fgrav — F cen t r if U g a i H- F e i the relation 

(15) U P {1 - q 2 ) = n 2 P . 

Thus one can infer the tangential derivative at the disk for given Cl and constant Q and then 
solve the boundary value problem for the Laplace equation as above. The electric and the 
gravitational potential are in this example just proportional. Similarly one can prescribe ad 
and obtain f l from {H}. Note that Q 2 has to be smaller than 1. For Q 2 = 1, the angular 
velocity in the dust disk vanishes. 

3. Einstein-Maxwell equations and group structure 

In this section we study the Einstein-Maxwell equations in the presence of one Killing 
vector. We explore the group structure of the equations and give the Harrison transforma¬ 
tion which generates electro-vacuum solutions from pure vacuum solutions. The solutions 
contain an additional real parameter related to the total charge. General properties of the 
transformed spacetimes as the asymptotics are discussed. In the stationary axisymmetric 
case, complete integrability of the equations is established. 

3.1. Maxwell equations. It is instructive to consider first the Maxwell equations in the 
absence of gravity, i.e. on a flat background. In standard notation, the equations for the 
electric and respectively magnetic fields E and B read 

(16) divE = 0, rotE + = 0, divB = 0, rotB — E >t = 0. 

Since divB = 0, one can define (up to gauge freedoms) a vector potential via B = rotA. 
It is convenient for the relativistic treatment in the following sections to introduce four¬ 
dimensional notation. We use the convention that greek indices take the values 0,1, 2,3 
and latin indices the values 1,2, 3. A four-dimensional vector potential A p = ( A , A a ) is 
introduced which is related to the tensor F pi/ of the electromagnetic fields via 

(17) F pv = A^ v - A Vift . 

In vacuum, the Maxwell equations read (indices are raised and lowered with the Minkowski 
metric) 

(18) F lu /“ = 0, *F^' V = ^e pva0 F a0 ’ v = 0. 

They can be obtained from the action 

(19) S= \j dxiF ^ F ^' 

Obviously the equations {24} are invariant under the discrete exchange F —> *F. In 
addition they are invariant under the continuous rotations 

(20) F + i*F->e ie (F + i*F), 9 e R, 

or in terms of the fields, E + zB —> e“(E + zB). This is the well known 17(1) symmetry 
of the Maxwell equations in vacuum. It can be used to ‘generate’ solutions from known 
ones: if for instance a solution with vanishing magnetic field is given, the f7(l)-symmetry 
can be applied to generate from the given electric field a solution with non-vanishing mag¬ 
netic field. The transformed fields contain a new real parameter. In case there is an electric 
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monopole moment, there will be a magnetic monopole moment in the transformed solu¬ 
tion. This often limits the physical relevance of the so generated solutions. 

In the stationary case, the potential /1 ; , can be chosen to be independent of the time 
coordinate t. Since rotB = 0, we can define a scalar potential via B = grad B, the well- 
known magnetic potential for stationary fields. The equation divB = 0 implies A B = 0. 
In four-dimensional language which will be needed in the Einstein-Maxwell case, this 
construction works as follows: Since F ab is a three-dimensional antisymmetric tensor, it 
can be dualized to an axial vector by contraction with the totally antisymmetric e-tensor. 
The Maxwell equations d F ah j, = 0 imply that this vector must be a gradient of some 
potential B. We can thus define the potential B via 

(21) B )C = ^abc^a.b- 

The potentials can be combined to the complex potential <I> = A + iB. In this case the 
action is just given by 

(22) S =^Jdx 3 |V$| 2 , 

and the Maxwell equations read 

(23) A$ = 0. 

In the stationary case, the Maxwell equations are thus equivalent to the Laplace equation 
for a single complex potential $. 

3.2. Einstein-Maxwell equations. In the Einstein-Maxwell case, the Maxwell equations 
have the same form as in d- only the partial derivatives have to be replaced by covariant 
derivatives since the spacetime is no longer flat, 

(24) F llv ,v = 0, *F^ V = 0. 

The tracefree energy-momentum tensor of the electromagnetic field is given by 

(25) r M „ = F m F v a - 
The Einstein equations have the form 

(26) R, lv - ^g^R = T I1V , 

where is the metric of the spacetime, R ln/ the Ricci tensor and R the Ricci scalar. With 
d we get for HD 

(27) R„ v = F^F* - iff ^F kX F kX . 

Equations d241i and d form the Einstein-Maxwell equations. They can be derived from 
the action 

(28) S = i / dx A ^g (r - ^F^F^ , g = det(g^). 

Since the Maxwell fields only enter the Einstein equations via the energy-momentum ten¬ 
sor, the C7(l) symmetry (1201 carries over to the Einstein-Maxwell case and can be used to 
generate solutions as before. 

In general one would expect that the above 17(1) invariance of the Einstein-Maxwell 
equations is the only symmetry of the equations even in the presence of Killing symmetries 
in the spacetime. However, it turns out that a much bigger symmetry group exists already 
for a single Killing vector. It is convenient to use a projection formalism which goes back 
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to Ehlers |i29|, see also 130| and 1381 for additional references. The case considered here 
can be viewed as a special case of the standard dimensional reduction of Kaluza-Klein 11381 
and supergravity theories 1391 . In this formalism, the metric is written in the form 

(29) ds 2 = —f(dt + k a dx a )(dt + k b dx b ) + jh. ab dx a dx b . 

We are considering here for convenience the stationary case, i.e. we use coordinates in 
which the timelike Killing vector is given by 9 t , all potentials are independent of t. How¬ 
ever, the results hold with minor changes for general Killing vectors. Coordinates are intro¬ 
duced for convenience, the dimensional reduction can also be carried out in a coordinate- 
independent way. We assume that the Killing vector is timelike throughout the spacetime, 
i.e. that its norm / does not vanish. 

The vector potential is decomposed as the metric into pieces parallel and orthogonal to 
the Killing vector, A M = (A, A m + k m A). The Lagrangian of (1281 can then be written in 
the form 

(30) 

C = ^Vh[n- ^h ab lal b + t-K,ab1C ab + - f h ab A, a A. b - j -(T ab + AYi ab ){T ab + AK, ab )^j 

where £ is a three-dimensional Lagrangian density, where T ab = A n j, — A b , a , and where 
1C a b = k a , b — fcb, a . All indices are raised and lowered with h ab . Note that the tensor JC ab 
vanishes only if the Killing vector is hypersurface orthogonal in which case the spacetime 
is static. 

The first part of the Maxwell equations (1241 can be written in the form 

( 31 ) -^(V=gFn, v = 0 

which implies (\/h.F ab //) tb = 0. With this relation or by varying Col with respect to A a , 
we obtain 

(32) (Vhf(F ab + AJC ab )) }b = 0. 

We can define as before the potential B via 

(33) B tC = -le cab Vhf(F ab + AJC ab ). 

Again A and B can be combined to the complex electromagnetic potential «1> = A + iB. 

Similarly we get by varying css with respect to k a 

(34) (Vft (^-KL ab - Af[T ab + AK, ab )^ =0. 

This can be dualized as above by introducing the so-called twist potential b via 

f 2 

(35) h, c = ^cab^fk—KA b + BA C — AB C . 

The potentials / and b can be combined to the complex Ernst potential, 

(36) £ = f - $!> + ib. 


The scalars b and B replace the vectors k a and A„. The corresponding three-dimensional 
Lagrangian reads with w a = 6 >a — 2 BA a + 2 AB a 

\[h 


(37) 


£ = 


n - h 


ab 


2 >b 


W a W b ) - j(A< 


Ah + B a B i 
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The line element 

(38) ds 2 = ^-((df) 2 + (db + 2BdA-2AdB) 2 - h(dA) 2 + (dB) 2 )) 

2 J- J 

describes the invariant metric of the Riemannian symmetric space S = SU (2, 1)/S[U(1, 1) x 
E/(l)] in some coordinates. The stationary Einstein-Maxwell equations can thus be in¬ 
terpreted as three-dimensional gravity coupled to some matter model. The ‘matter’ is a 
SU{2, l)/5[t/(l, 1) x E/(l)] nonlinear sigma model 140114II . Note that sigma models are 
related to harmonic maps P42l. The space S can be parametrized by trigonal 3x3 matrices 
V, 

/ v7 o o \ 

(39) V = iV 2$ _1 ° 

V (& + *l$| 2 )A/7 (V2$)/77 i/v7 / 

The matrix V satisfies 

/ 0 0 i \ 

(40) V^rjV = 77, 77 = 0 1 0 , 

\ -i 0 0 / 

i.e. it is unitary with respect to the metric 77 of SU( 2,1). The action of Q £ SU(2, 1) on 
V is 

(41) V —> H(V,Q)VQ~ 1 , H(V,g)eS[U(l,l)xU(l)], 

where H restores the triangular gauge of V. To obtain a gauge invariant parametrization, 
one introduces 

(42) X := EVtSV, S = diag(l, -1,1), 
on which the action of Q £ SU (2,1) is given by 

(43) x —► ZiQ-^SxS- 1 . 

We have 

(44) 

( f ~ 2|$| 2 + (b 2 + | < f>| 4 )// V 2 Hb - *|4 >| 2 + */)// (b - m 2 )/f \ 

X= -\/2$(&-H|$| 2 -i/)// 1 — 2|4>| 2 // -y/2 $// . 

V (b + im/f y/2$/f 1 // / 

The SU (2,1) symmetry can be used to generate solutions by the action of an element Q. 
We list the infinitesimal transformations and their consequences: 

/ 0 0 0 \ 

\9i 0 0 

\ 02 03 0 J 

lead to gauge transformations which add physically irrelevant constants to A£ and 

/ 0 0 9 \ 

0 0 0 
\ 0 0 0 / 

is an Ehlers transformation ia which changes / —> b, i.e. which generates station¬ 
ary from static solutions (if the ADM-mass of the spacetime is non-zero, the transformed 
solution will have a Newman-Unti-Tambourini (NUT) parameter which corresponds to a 
magnetic monopole and which is believed to be unphysical) 
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iO 0 0 

0 -2 iO 0 

oo ie 

is an electromagnetic duality transformation, i.e. the U( 1) symmetry of the sourcefree 
Maxwell equations, 



is a scale transformation, /, b, $ —> e 6 /, e B b, e e / 2 <f>, and 
0 -id 0 
0 0 0 
0 0 0 

is a Harrison transformation ED which changes / —> $, i.e. generates solutions with 
electromagnetic fields from pure vacuum solutions. 


3.3. Harrison transformations. Harrison transformations offer the possibility to gen¬ 
erate solutions with charge from pure vacuum solution. This leads to solutions to the 
Einstein-Maxwell equations containing one additional constant parameter which is related 
to the charge. For physical reasons we are interested in solutions which are equatorially 
symmetric and asymptotically flat, i.e. / —> 1, d> —> 0 and 6^0 for |£| —» oo. Asymptotic 
flatness is just the mathematical formulation of the physical concept of an isolated matter 
distribution, e.g. a galaxy. Solutions with equatorial symmetry, i.e. a class of solutions 
where the metric functions have a reflection symmetry at the equatorial plane ( = 0, are 
of special physical interest. In a Newtonian setting it can be proven that perfect fluids in 
thermodynamical equilibrium lead to equatorially symmetric situations, and the same is 
assumed to hold in a general relativistic context. A consequence of this condition is that 
NUT-parameters are ruled out. 

We assume that the pure vacuum solutions which we want to submit to a Harrison 
transformation satisfy these conditions. To ensure that the transformed solutions have the 
same asymptotic behavior, one has to use a scale transformation (/ —> 1) together with a 
transformation which changes $ and b by some constant ($, 6 —> 0). By exponentiating the 
matrices of the SU(2 ,1) transformations, we thus consider a transformation of the form 


(1 

id i 

—z0i0i/2 \ / 


0 

0 \ 

(45) 0=0 

1 

-0i ] diag(e~ 92 , 1, e® 2 ) | 

-03 

1 

0 

V o 

0 

1 J \ 

— 04 + 0305/2 

— 05 

1/ 


Since the asymptotic conditions imply that \ and x' are the unit matrix at infinity, the 
matrix Q must satisfy the condition EC/'E = Q~ 1 . This leads with 1451 to 

(46) e e *= t ^ 03 = *0i, 04 = 0, 05 = 01- 

The matrix Q thus takes the form 

1 i0i_ -id i<9i/2 \ 

— i6\ 1 + didi/2 — 6\ I . 

ie i0”i/2 -01 1 J 


(47) 




1 - 0 i 0 i /2 


If we transform an Ernst potential of a pure vacuum solution (<f> = 0), we end up with 

f 48 t 6 ^~t (/ 2 + &2 )/ 2 ~ + W 2 )/ + 1 - Mi - 0 i 0 ~i/ 2 ) 

1 J V2 (1 - 0 ! 0 " l // 2)2 + ( 0 1 0 ' 1 / 2)202 
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We are interested in transformations which preserve the equatorial symmetry, i.e. /(—C) = 
/(£), 6(— C) = — &(C) and $(— C) = 3>(C)- This implies for d4^t that 9\ must be real 
which rules out magnetic monopoles. We put q = 0\ /\/2 and sum up the results for the 
transformed potentials: 


(49) 

(50) 

(51) 


(i-g 2 ) 2 / 

(1 — q 2 f ) 2 + q 4 b 2 ’ 

(1 ~q 4 )b 

(1 — q 2 f) 2 + q 4 b 2 ’ 

(1 - f)(l - q 2 f) + q 2 b 2 + ib(l - q 2 ) 
q (1 — q 2 f) 2 + q 4 b 2 


The real parameter q has to be in the region () < \q\ < 1, for q > 1 the transformed 
spacetime would have a negative mass if the original mass was positive. The value q = 0 
corresponds to the untransformed solution. The above formulas imply that the functions 
/', b' and +' are analytic where the original functions are analytic. 

A well-known example is the Harrison transformation of the Schwarzschild solution 
which leads to the Reissner-Nordstrom solution. In the Ernst picture the Schwarzschild 
solution reads in cylindrical Weyl coordinates with r± = ± to) 2 + p 2 


(52) 


/ = 


r_|_ + — 2 to 

r_|_ + r_ + 2 to ’ 


6 = 0 , 


where the horizon (/ = 0) is located on the axis between — m and m. For the transformed 
solution we get with <El) 

(53) 

, = (r + + r_) 2 — 4m 2 , = 2Q , = 1 + q 2 = 2mq 

(r + + r_ + 2m') 2 ’ r+ + r_ + 2m' ’ 1 — q 2 ’ 1 — q 2 ’ 


and b' = 0 which is the Reissner-Nordstrom solution. This is a static spacetime with mass 
m' and charge Q subject to the relation m' 2 — Q 2 = m 2 . Both m! and Q diverge for 
q —> 1. The extreme Reissner-Nordstrom solution with m! = Q is only possible in the 
limit to —> 0, | q \ —> 1. The horizon of the solution is again located on the axis between 
±m which illustrates that the horizon degenerates in the extreme case. 


3.4. Asymptotic behavior of the Harrison transformed solutions. We assume that the 
asymptotic behavior of the original solution, which can be read off on the axis, is of the 
form / = 1 — 2M/|C|, 6 = —2 J/C, 2 and 4> = Q/\C,\ — IJm/C, 2 pl us terms of lower 
order in 1/|C| where M is the Arnowitt-Deser-Misner mass, J the angular momentum, Q 
the electric charge and Jm the magnetic moment. The same will hold for the Harrison 
transformed potentials. We find |2(J| 


(54) 


and 


(55) 


M' = M 


1 + q 2 


2 q 


1 -q 2 1 -q- 


;Q, J' = J 


1+q 2 2 q 

1 -q 2 1 -q- 


■Jm, 


Q' = Q 


1 + q 2 2 q 

1-q 2 1 -q- 


■M , — Jm 


1 + q 2 2 q 

1-q 2 1-q 2 


:J. 


It is interesting to note that the quantities M 2 — Q 2 and ,/ 2 — J 2 M are invariants of the 
transformation. They are related to the Casimir operator of the SU( 2, l)-group. If the 
original solution was uncharged, the extreme relation M' = +(}' is only possible in the 
limit M -> 0. 



















CHARGED DUST DISKS 


13 


A consequence of the relations is the presence of a non-vanishing charge if the 
ADM mass of the original solution is non-zero whereas the charge is. Since charges 
normally compensate in astrophysical settings, this limits the astrophysical relevance of 
Harrison transformed solutions. 

A further invariant is the combination JmM — JQ which is of importance in relation 
to the gyromagnetic ratio 


(56) 


2 MJm 
9m = -jq- 


Relation implies that g' M is equal to 2 if Q = Jm = 0 and q / 0. Thus all solu¬ 
tions which can be generated via a Harrison transformation from solutions with vanishing 
electromagnetic fields as the Kerr-Newman family from Kerr have a gyromagnetic ratio of 
2 . Due to the invariance of JmM — JQ under Harrison transformations, a gyromagnetic 
ratio of 2 is not changed under the transformation. 

Whether this property is an indication of a deep relation between relativistic quantum 
mechanics and general relativity as claimed in |27 [ 28 | is an open question. Here it is 
just related to an invariant of the Harrison transformation, a subgroup of SU( 2,1). Since 
most of the known solutions to the Einstein-Maxwell equations can be generated via a 
Harrison transformation from solutions to the pure vacuum equations as the Kerr-Newman 
family from Kerr, a gyromagnetic ratio of 2 is well known from exact solutions. Numerical 
calculations of charged neutron stars E3 indicate, however, that values well below 2 are 
to be expected in astrophysically realistic situations. 


3.5. The stationary axisymmetric case. In the astrophysically important stationary ax- 
isymmetric case, the symmetry group of the equations increases again, this time to the 
infinite dimensional Geroch group 1 44. 45 1 . This means that the equations are completely 
integrable, where the notion of integrability is to be understood in a Hamiltonian sense: 
the equations have the same number of conserved quantities as degrees of freedom. For 
completely integrable systems, this number tends in a countable way to infinity. The in¬ 
finite dimensional symmetry group shows up in treating the differential equation under 
consideration as the integrability condition of an overdetermined linear differential system 
for a matrix-valued function. The system contains an additional parameter, the so-called 
spectral parameter which reflects the infinite dimensional symmetry group. 

In the presence of a second Killing vector, the metric J29b can be further specialized. 
The axial Killing vector d 0 commutes with the timelike Killing vector Of . The metric h a b 
of J29b can be chosen to be diagonal, h a b = diag(e 2fe , e 2fc , p 2 ), and k a can be brought into 
the form k a = (0, 0, a) which leads to the Weyl-Lewis-Papapetrou metric, see e.g. |46]| 

(57) ds 2 = — e 2U (dt + ad(p) 2 + e~ 2U (e 2fc (d p 2 + d£ 2 ) + p 2 d(/) 2 ) , 

where / = e 2U . In this case the Einstein-Maxwell equations reduce to 

fA£ = (V£ + 2<I>V<I>)V£, 

(58) /A<1> = (V£ + 2i>V$)V$, 

where A and V are the standard differential operators in cylindrical coordinates, and where 
the potentials £ and $ are independent of (j >. The first equation generalizes the Newto¬ 
nian equation A U = 0 to the Einstein-Maxwell case, the second equation generalizes the 
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Maxwell equation (1231 . The duality relations d33l and d35l read 


(59) 

(60) 


—a,A t £), 

H = (i(^)£ + ~ $$«)» 


which implies that a and Ag follow from 8 and <1'. We choose a gauge where A\ = A 2 = 0. 
The equations for TZ a i, of d37) are equivalent to 


(61) 


k = 


Z -1 ( i 




Thus the complete metric and the electromagnetic potential can be obtained from given 
potentials £ and $ via quadratures. 

The system d58> was shown to be completely integrable in |24j. In the form 1471 . the 
associated linear differential system for a 3 x 3 matrix-valued function T reads (we use the 
complex coordinate £ = ( — ip) 


(62^ftf _1 


Pi 0 Adi \ K p ( 0 Pi 0 \ 

0 Ci 0 +-- Ci 0 -Adi , 

-A/i 0 i(Ci+Pi) ) V 0 -A/i 0 ) 

P 2 0 Ad 2 \ K -8 ( 0 V 2 0 \ 

0 C 2 0 +- i C 2 0 -Ad 2 , 

-U 2 o i(C 2 + P 2 ) ) V 0 -U 2 o ) 


where T depends on the spectral parameter K which varies on the Riemann surface £ 
of genus zero given by the relation //q( K ) = (K — 8){K — £). This is a first hint on 
the relevance of Riemann surfaces in this context. Notice the special feature of the Ernst 
equation that the branch points £ £ depend on the spacetime coordinates. We denote a 
point on £ with the projection K in the complex plane by P = (K , ±//,o(A')) = A . On 
£ there is an involution <j that interchanges the sheets, i.e. with P = (A', /ro(A')) we have 
aP = P a = (A', —fia(K)). We use the notation ocA for the infinite points on different 
sheets of the curve £, namely n/K 9+l —> ±1 as K —> oA. 

The expressions for Cj, P; and Ad , (i = 1,2) follow from the condition 


(63) 


( £ + 2d>$ 

1 72i<r ^ 

1 ( 1 

0 

V2i4> 

^(ocA, z.z)= ( 8 

-1 -v^i® 

= 0 

1 

0 

\ -2i$A 

0 V2e u ) 

V o 

0 

V2e u 


Thus we have det 'T(ocA) = — 2\^2e 3U . The reality conditions read 


8 1 

8 -1 

-V2id> 0 


0 

-V2id- 

1 


(64) C 2 = Pi, Ci=P 2 , Af 2 = -Mi, Ad 2 = —U 2 . 
The inverse matrix takes the form 

/ 1 1 0 \ 

(65) Tr-1 = £ + -£_ -2 i<f>e u , 

\ y/2 x/2 V2e u ) 
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which implies via (1621 

Ci 

= (% + 2#^)/(2/), 


Pi 

= (C 5 + 2$$ ? )/(2/), 


Adi 

= i<p£e~ u , 

( 66 ) 

A/i 

= i^e~ ^ . 


The system © can be used to generate solutions to the Ernst equations. To this end 
one investigates the singularity structure of the matrices and with respect 

to the spectral parameter and infers a set of conditions for the matrix T which satisfies 
the linear system d62t . These conditions can be summarized in the following theorem (see 

My 


Theorem 3.1. The matrix T is at least twice differentiable with respect to £ and £. T (P) 
is holomorphic and invertible at the branch points !; and £ such that the logarithmic deriv¬ 
ative has a pole at £ and has a pole at £. 

II. All singularities of'll on C are such that the logarithmic derivatives 1 and T 

are holomorphic there. 

III. The matrix function 'T is subject to the reduction condition 

(67) eT>(P CT ) = *(P) 7 (P), 

where e = diag( 1 , — 1 , 1 ), and where 7 is an invertible matrix independent of £, £. 

IV. At oo + , the matrix function is given by \63\ . 

A proof of this theorem may be obtained by comparing the matrix T with the linear 
system d62t . 

Proof. Because of I, T and T 1 can be expanded in a series in the local parameters tv = 
s/K — £ and T| = \JK — £ in a neighborhood of P = £ and P = (, f tj respectively at 
all points £, £ which are not singularities of d/. This implies that d^d / -1 = ao/f + + 

ct 2 t + .... We recognize that, because of I and II, '^^'h~ 1 — ao/t is a holomorphic function. 
The normalization condition IV implies that this quantity is bounded at infinity. According 
to Liouville’s theorem, it is a constant. Since d/, T “ 1 and d 7 are single valued functions 

on £, they must be functions of K and p^. Therefore we have d/jd ) -1 = + 

(3\. The matrix If must be independent of K and po since d/jd / -1 must have the same 
number of zeros and poles on C. The structure of the matrices If and /3i follows from III. 
From the normalization condition IV, it follows that has the structure of ©■ The 

corresponding equation for 'l7 'l ; 1 can be obtained in the same way. □ 

The choice of the matrix 7 in d8 II corresponds to a gauge freedom which is however not 
completely fixed. A matrix C(K) with the property [C, 7 ] = 0 with C(oo) = 1 which acts 
on T in the form d/ — >■ '1>C( K) leads again to a solution of (1621 for the same potentials £ 
and d>. 

We note that the metric function a can be directly obtained from the matrix T without 
integrating ©: We denote by DpF(P) the coefficient of the linear term in the expansion 
of the function F in the local parameter near P. Using the identity (d ) -3 I)^+ = 

d l ~ 1 D 00 + (dyd'~ 1 )d', one finds with ( 1621 to (1661 

( 68 ) ((’T- 1 P 0O+ T')b 12 = -^(Cr - Vf) = -^(i(9£) £ + d>dy - $$*). 
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Equation (1601 leads to the expression 

(69) (a — ao)e 2C/ = *Z3 00 +('I'i2 + (t , 22 ), 

which relates the metric function a directly to the potential >1-'; the constant a o is determined 
by the condition that a has to vanish on the regular part of the axis. 

To generate solutions to the Ernst equation one has to construct matrices according to the 
above theorem. The theorem ensures that the so constructed solutions to the Ernst equation 
are analytic at all points where the conditions I-IV are satisfied. It also makes clear where 
singularities of the spacetime can be found: if the moving branch points £, £ coincide with 
singularities of \k, condition I is no longer fulfilled and the theorem does not hold. These 
points are possible singularities, but they could be regular if the singularities of >h are for 
instance pure gauge. On the other hand, if the singularities of the spacetime are prescribed 
— for electro-vacuum the boundary of some matter source constitutes a singularity — one 
has to construct a matrix with the corresponding singularity. This means one has to solve a 
so-called Riemann-Hilbert problem, i.e. to find a function with prescribed singularities on 
C. 

This method to construct solutions to integrable equations is also known as the inverse 
scattering method |48l . Explicit solutions are in general only known for scalar Riemann- 
Hilbert problems which lead in the case of the Ernst equation to static solutions for the 
pure vacuum, see 114111511761 . In this case the Ernst equation reduces to the axisymmet- 
ric Laplace equation, the Euler-Darboux equation. In the matrix case, the solution to a 
Riemann-Hilbert problem is equivalent to an integral equation, see EHEma for the 
Ernst equation. A simple special case are so-called ‘soliton’-solutions or Backhand trans¬ 
formations where the matrix ok has only poles and zeros of the determinant as singularities. 
Since the Ernst equation is an elliptic equation, it has obviously no solutions describing 
physical solitons, but it has solutions with the same mathematical properties. The most 
prominent representant of this class for the Ernst equation is the Kerr-Newman family of 
charged rotating black holes. 

In na it was shown for the pure vacuum case that the Riemann-Hilbert problem can 
be solved explicitly for a large class of problems by exploiting the gauge freedom of the 
matrix ok. This leads to solutions which are defined on certain Riemann surfaces. It is well 
known that large classes of solutions to non-linear integrable equations can be constructed 
via methods from algebraic geometry. For evolution equations like Korteweg-de Vries 
and Sine-Gordon, see e.g. 1501 . the corresponding solutions are periodic or quasi-periodic. 
Algebro-geometric solutions to the Ernst equation do not show any periodicity as will be 
discussed in the following section. A way to construct such solutions is via the so-called 
monodromy matrix [51 501. For a linear system of the form 

(70) ^ = kVk, >k 5 - = V\k 

as m, the monodromy matrix L can be defined as the solution of the linear differential 
system 

(71) L t = [W,L\, L|=[y,L]. 

It follows from equation <EB that the characteristic polynomial 

(72) Q(/i, K) = det(L - jll) 

is independent of the physical coordinates (the coefficients of the polynomial are ‘integrals 
of motion’). The equation Q(/t, K) = 0 is then the equation of a plane algebraic curve. 
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Since L is a 3 x 3-matrix, relation M2\ is cubic in /'/ which can be always brought into 
normal form by a redefinition of fi : 

(73) fi 3 +V(K)ji + Q(K) = 0. 

The functions V and Q are analytic in K. Equation <1731 defines a three-sheeted Riemann 
surface which will in general have infinite genus. For polynomial V and Q, the surface 
will be compact and will have finite genus. For a given surface the solutions to the Ernst 
equations can be given in terms of the theta functions on this surface which was first done 
in fj~4l for a special case. 

Since the theory of these surfaces is not as well understood as the theory of hyperelliptic 
surfaces which occur in the pure vacuum case, we will study in the next section Harrison 
transformed hyperelliptic solutions. In section 6, we will come back to solutions to the 
Ernst equations on three-sheeted surfaces. 

3.6. Conformastationary and magnetostatic solutions. There are two special cases of 
the above equations for which disk sources have been considered: Conformastationary 
metrics |52 1531 are included in the above formalism as metrics with a flat metric h a b, 
there is no axial symmetry required. Writing / = (VV)~ l , the function V is a complex 
solution to the three-dimensional Laplace equation. For asymptotically flat solutions V has 
to tend to 1 at infinity. The electromagnetic fields are given by 

(74) E + *H = grad , D + *B = |V|(E + £H) + iTx (E + zH), 
where T is a solution to 

(75) rotT = i(VgradV — LgradV). 

Static spacetimes with a pure magnetic field can be mapped to the case of stationary 
axisymmetric vacuum spacetimes, see 1461 [Til. Writing A^ = (0,0, 0, A(p, £)), the field 
equations reduce to the Ernst equation which is discussed in the next section. Thus all 
solutions to the vacuum Ernst equation can be interpreted as magnetostatic solutions by 
putting / — » l/(9?£ ) 2 and A —> — V2?s£. 

4. Stationary axisymmetric Einstein equations, theta functional 

SOLUTIONS AND COUNTER-ROTATING DUST DISKS UNDER HARRISON 
TRANSFORMATIONS 

In this section we will summarize results on the Ernst equation in the pure vacuum 
case, a class of solutions on hyperelliptic Riemann surfaces m and a member of this 
class describing counter-rotating dust disks EU which were discussed in E3 and ED- 
We study the action of a Harrison transformation on stationary axisymmetric solutions and 
discuss the transformed counter-rotating dust disk as an example. 

4.1. The stationary axisymmetric vacuum. In the stationary axisymmetric vacuum, the 
Einstein-Maxwell equations of the previous section hold with $ = 0. Since the simplifica¬ 
tion is considerable, we list the relevant equations below. The Ernst potential £ = / + ib 
is subject to the Ernst equation [55 ] 

<76) - 2(^T) <£ f - £i) = 
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If the Ernst potential is real, equation (1761 is equivalent to AU = 0. The corresponding 
spacetime is static and belongs to the so-called Weyl class. For a given Ernst potential the 
metric follows from 


(77) 

(78) 


a t 


k 


2 P 


(£ - kt 

(£ + £) 2 ’ 

kk 


(e-o 


(£ + £) 2 ' 


The stationary axisymmetric Einstein equations in vacuum were shown to be completely 
integrable in [56] and B71 . The associated linear system can be formulated for a 2 x 2- 
matrix ^ (see 1581 ) 




-l 


(79) 
where 

(80) 


'I'r'f 


-1 


M 0 

o a f 
a7 o 

o M 



K-Z ( 0 M 


K-£\ Af 0 

K-if 0 M 
K-i\ M 0 


M = 


k 


Af = 


£ + £’ £4 

Theorem 3.1 holds with the following changes: 

111’. T is subject to the reduction condition 

(81) 'f(P a ) = a 3 ^{P)a 2 , 

where <j 2 , 03 are Pauli matrices. 

IV’. The normalization and reality condition 


(82) 


Tqp = oo 1 ) = 


£ —* 
£ i 


The function £ in \82\ is then a solution to the Ernst equation dm 

We note that the choice of the normalization of T at infinity is different from the one 
in theorem 3.1. This form was chosen to implement the Harrison transformation as was 
done in the last section with the metric 77 in ( 1401 of SU{2, 1). The choice of the matrix 

0 -i 


02 = 


0 


in condition III is again a gauge condition, which does not fix the 


gauge completely, however. The remaining gauge freedom is here due to matrices C(K) = 
ni{I\ )1 +k 2 (K)ct 2 , where the n t do not depend on £, £ and obey the asymptotic conditions 
Ki(oo) = 1 and n 2 (oo) = 0. The matrices C act on in the form T —> ^C(K). Since it 
is a consequence of ( 1791 that det 'I' = F(I\)e 2U where F(K) is independent of £, £ we can 
use this gauge freedom to choose F[K) = 1. The linear system ( 1791 leads for the matrix 
X(P) = ’T _1 (oo _ )T'(P) to one of the linear systems used in 1451 . This parametrization, 
especially 


(83) 


T ^00 )'3/(oo + ) = 


2 ££ *(£ - £) 
*(£-£) 2 


reveals that the Ernst equation is an SL(2,S.)/SO(2) sigma model. For more details on 
the group aspect see |45l and the discussion in the previous section. 
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For given 'k(/i), one can again directly determine the metric function a without having 
to integrate relation (I77> . With < 1791 one finds for the matrix S = \I/~ 1 (oo + )D 00 +'I , (oo + ) 


,84) S f := (t-'foo +)D^) t = 

This implies with G3 

(85) a — oq = — 2£>i2 = — 2(\D 1 (oo + )I7 00 + v I , )i2, 


«£-£)* \ 
-(ssh )' 


where ao is a constant which is fixed by the condition that a vanishes on the regular part 
of the axis. The matrix S will be needed to calculate the metric function a after a Harrison 
transformation. 

The monodromy matrix L defined in {H} corresponding to G3 is now also a 2 x 2- 
matrix. Since this matrix can be chosen without loss of generality to be tracefree, the 
characteristic equation <EJ takes the form 

(86) fi 2 = V(K), 

where V is an analytic function in K. For polynomial V = YFi=i (K ~ F) )( K — Fi), the 
two-sheeted Riemann surface defined by 1861 is compact. Since the spectral parameter K 
varies on the two-sheeted surface C, equation <1861 defines a two-sheeted branched cover of 
C. In other words the so defined Riemann surface £ is a four-sheeted cover of the complex 
plane. The structure of this surface is shown in the Hurwitz diagram Fig. [2 















E 

;f e 

T •" e 

9 * 

~'(T 

9 











C £ El F 1 ■■■ Eg Fg 


Figure 1. The Hurwitz diagram of £ shows the Riemann surface as 
seen from the side. 

There is an automorphism er of £ inherited from £ which ensures E?, Ei and F) 7 , 
Fi have the same projection in the complex plane. The orbit space £h = £/cr is then, 
see 1501 . again a Riemann surface given by 

g 

(87) p 2 = (K - 0(K - 0 - Ei)(K - Fi). 

The fixed points £, £ of the involution a are additional branch points of the Riemann sur¬ 
face. The points Ei — ai — id, have to be constant with respect to the physical coordinates. 
They are subject to the reality condition E, = F t or Ej,F % e I for i = 1,,g. A sur¬ 
face of the form <|87j is called hyperelliptic, since the square root of a polynomial can be 
considered as the straight forward generalization of elliptic surfaces. Thus it is possible to 
construct components of the matrix T on £ // which makes it possible to use the powerful 
calculus of hyperelliptic Riemann surfaces. 
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4.2. Solutions on hyperelliptic surfaces. Solutions to the Ernst equation on hyperelliptic 
surfaces were given by Korotkin in d. and in the gauge (1821 in 1581 . We summarize 
basic facts of hyperelliptic surfaces (see e.g. |50) and (60) to ED) to be able to present 
these solutions. 

A Riemann surface has a topological invariant, the genus, which loosely speaking gives 
the number of holes in the surface. A surface of genus g > 1 is topologically just a sphere 
with g handles. For genus 0, it is the Riemann sphere, in the elliptic case g = la torus. 
The hyperelliptic surface Ch defined by has genus g. We order the branch points with 
Im Ei < 0 in a way that R/V| < KP 2 < ... < $tE g and assume for simplicity that the real 
parts of the Ei are all different, and that there are no real branch points. On this surface we 
introduce a canonical basis of cycles which are non-homologous to zero, i.e. which cannot 
be contracted to a point. This basis consists of 2 g cycles a;, bi, i = 1,..., g which do not 
intersect except for a,;, bi with the same index. For the surface of genus 2 we will consider 
in detail later, we use the cut-system in Fig. [2] The surface resulting from cutting Ch 



Figure 2. Canonical cycles (Pq = 0. 


along these cycles, the fundamental polygon, is simply connected. 

On a surface of genus g , there are g independent holomorphic one-forms which are also 
called differentials of the first kind. These one-forms can be locally written as F(K)dI\ 
where F(K) is a holomorphic function. Their integrals are holomorphic functions. For 
general Riemann surfaces, the determination of the holomorphic differentials is a non¬ 
trivial problem. One important simplification in the theory of hyperelliptic surfaces is that 
these differentials are explicitly known. A basis is provided by 

,'cLK KdK K 9 ~ 1 dK' 

( 88 ) dv k = -, 


MM M 

The holomorphic differentials du>k are normalized by the condition on the a-periods 


(89) 


dui k — ‘I'lriSik. 


The matrix of 6 -periods is given by B^ = j h dujk■ The matrix B is a so-called Riemann 
matrix, i.e. it is symmetric and has a negative definite real part. The Abel map ui : Ch —■* 
Jac (Ch) with base point E\ is defined as w(P) = diOk, where Jac (Ch) is the Jacobian 
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of Ch , factorized with respect to the lattice of a- and 6 -periods. The theta function with 
characteristics corresponding to the curve Ch is given by 

(90) e pq (x|B) = exp | ^(B(p + n), (p + n)) + (p + n, + x)| , 

n£Z9 ' ' 


where x £ C 9 is the argument and p,q e C s are the characteristics. We will mainly 
consider half-integer characteristics in the following. A half-integer characteristic is called 
odd if 4(p,q) is odd, and even if this inner product is even. The theta function with 
characteristics is, up to an exponential factor, equivalent to the theta function with zero 
characteristic (the Riemann theta function is denoted with 0) and shifted argument, 

(91) 0 pq (x|B) = 0(x + Bp + 27Tiq)exp|i(Bp,p) + (p, 27riq + x)| . 

The theta function has the periodicity properties 

(92) 0 pq (z + 27 Tiej) = e 2 "^0 pq (z), 0 pq (z + Be,-) = 0 p q(z), 


where e, is the g-dimensional vector consisting of zeros except for a 1 in jth position. 

We denote by du/pQ a differential of the third kind, i.e., a one-form which is holo- 
morphic except for two poles in P,Q £ Ch with residues +1 and —1 respectively. This 
singularity structure characterizes the differentials only up to holomorphic differentials. 
They can be uniquely determined by the normalization condition that all a-periods vanish. 
The differential duj^ c + oc - is given up to holomorphic differentials by — K 9 dK) (i. 

In W7II18I a physically interesting subclass of Korotkin’s solution was identified which 
can be written in the form 


(93) 


0pq(u;(oo + ) + U ) J 

0pq(u;(oo-) + u) 


where the characteristic is subject to the reality condition Bp + q £ zM, where u = (u k ) £ 
C 3 and where 


(94) I = -T f In G(K) du> 00 + 00 - (K), u k =f lnG(K)dw k . 

Z7TI J p Z7TI J p 

T is a piece-wise smooth contour on Ch and G(K) is a non-zero Holder-continuous func¬ 
tion on T. The contour T and the function G have to satisfy the reality conditions that 
with K £ T also K £ T and G(K) = G(K); both are independent of the physical co¬ 
ordinates. In the case of disks of radius 1 in which we are interested here, the contour T 
is the covering of the imaginary axis in the +-sheet of Ch between —i and i. In ITSl it 
was shown that solutions of the above form on a Riemann surface of even genus g = 2s 
given by g, 2 = ( K — £)(K — £) rii=i(A 2 — — E 2 ) with a function G subject to 

G(—K) = G(K) lead to an equatorially symmetric Ernst potential, £(() = £(—£). 

Notice that these solutions depend only via the branch points of Ch on the physical 
coordinates. For 5 = 0 there are no theta functions and the potential (I93> takes the form 


(95) 


ln£ = 


27TZ 


In GdK 

V(K - C) 2 + p 2 ' 


As shown in section 2, this Ernst potential is equivalent to the Poisson integral with a dis¬ 
tributional density at the disk. The solutions can be considered as the static or Newtonian 
limit of the theta functional solutions 

In the theta functional solutions of evolution equations, the underlying Riemann surface 
is independent of the physical coordinates. The coordinates only occur in the argument of 
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the theta functions. Because of d92l the solutions are periodic or quasi-periodic. In the case 
of the Ernst equation, the Riemann surface itself is ‘dynamical’ since some of its branch 
points are parametrized by the physical coordinates. Thus the solution is given on a whole 
family of surfaces. The argument of the theta functions in ( 1931 has however 

no additional dependence on the physical coordinates. Here the modular dependence of 
the theta functions is important. Therefore the solutions show no periodicity and can be 
asymptotically flat. 

The limit £) —> F, in which the cut collapses corresponds to the solitonic limit for 
theta functional solutions to nonlinear evolution equations. The almost periodic solutions 
can be seen as an infinite train of solitons. In the ‘solitonic’ limit, only a finite number of 
solitons survive. In the case of the Ernst equation, the Kerr solution can be obtained as 
the corresponding limit of a genus 2 surface, see M- In this mathematical sense, black 
holes can be considered as solitons. The defining equation < 1871 for the Riemann surface 
illustrates the relation between solitons and theta functional solutions: the solitons have a 
different monodromy property which can be obtained in the limit of a degenerate surface 
(collapsing cut). Zeros of first order of V in <1861 correspond to theta functional solutions, 
zeros of order 2 to solitons. The 7V-solitons ea and the Backlund transformations are thus 
contained in Korotkin’s solutions as a limiting case. 

The importance of many exact solutions to Einstein’s equations (see |46|) is somewhat 
limited by the singularities they can have. In the stationary vacuum, the Lichnerowicz 
theorem 1631 states that a solution is either Minkowski spacetime or it has singularities. 
Non-trivial solutions must have sources which show in the form of singularities in matter- 
free settings. For general explicit solutions to the Einstein equations, it is however difficult 
to localize the singularities and to show that they can be replaced by matter sources. This 
is also in general true for Korotkin’s solutions. An important advantage of the subclass 
discussed here is that general statements can be made on the regularity of the solutions: 


Theorem 4.1. The Ernst potential (E3 with [pq] being a non-singular half-integer char¬ 
acteristic is analytic in the exterior of the disk £ = 0 , 0 < p < 1 where it is discontinuous 

iff 

(96) 0 pq (w(oo - ) + u) 0. 


For a proof see CD- Condition J96b reflects the fact that within general relativity ar¬ 
bitrary amounts of energy cannot be concentrated in a finite region of spacetime without 
forming a black hole or a singularity. This defines the range of the physical parameters 
where the solution is regular. A general method how to define this parameter range based 
on a study of the zeroes of theta functions was developed in l32l 

The Ernst potential J93l follows with (f82l> and J = 2^7 Jr In Gdupp* from a matrix + 
of the form 


(97) 

= e //2 


det(oo) 

det(A') 


e pq (u-MP)) j/2 

0 pq (u+u(oo _ ))_ 

©pq(u + +P)+++) J/2 

©pq (u+w(oo~)+u;(£)) 


_• Qpq (u+uj(P j )) j/2 

4 0 pq (u + w(c»-)) 

■ e pq (U + w(P g )+^(g)) -J/2 

©pq(u+w(00-)+^(C)) 


with (let (A) = 0 pq (u + u2(P))0 pq (u - u(P) + w(£)) + 0 pq (u - w(P))0 pq (u + 
w(P) +tu(0). 

To determine the metric function a via m, we have to calculate the matrix S which 
leads with d97l to 


<Sl 2 — 


2 if 


-D^+ In ■ 


0 pq (u + 0;(00 )) 
0 pq (u + w(oo-) +W0)’ 


(98) 
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and to determine the Harrison transformed function a in the next section. 


(99) 


8S_ e pq (u-!-^(00+)i 

2 */ °° + 1 0 pq (u + w(oc+) + a; (I)) 


Using a degenerated version of Fay’s trisecant identity ||62j (see Dl for the present case), 
we can write the above relations free of derivatives. 


P_ f Qpq(u)Qpq(u + 2tu(0O ) + ^(Q) _ 

2/ \U0 pq (u + w(oo-) + w(f))©pq(u + uj( 00 “)) 

p££ ( 6pq(u)Qpq(u + 2u;(oo + ) + cu(g)) 

2/ VL0 pq (u + w(oo+) + w(O)0 pq (u + w(oo+)) 

_ 0(w(oo _ ))0(w(oo _ ) + w(£)) 

" 0(O)©(w(O) • 

To construct the solution for the counter-rotating disks in we used an algebraic ap¬ 
proach which made it possible to establish algebraic relations between the metric functions 
at the disk. Let us recall that a divisor X on Ch is a formal symbol X = n\P\ +.. .+nuPk 
with Pi £ Ch and m € Z. The degree of a divisor is 2 j i=1 n i- The Riemann vector I\ r is 
defined by the condition that 0(u;(W) + Kr) = 0 if W is a divisor of degree g 1 or less. 
We use here and in the following the notation u>(W) = J Po dtu = J2i=i w ( Wi ). Note 
that the Riemann vector can be expressed through half-periods in the case of a hyperellip- 
tic surface. We define the divisor X = Yld -1 -T, as the solution of the Jacobi inversion 
problem (■£ = !,...,</) 




«Sl 2 

( 100 ) 

<^21 

where 


( 101 ) 



( 102 ) 


ui(X) — u)(D) = u, 


where the divisor D = (this corresponds to a choice of the characteristic in J93t > ). 

With the help of these divisors, we can write J931 > in the form 


r x 


(103) 


Inf = 


ID 


t 9 dr 

P-( T ) 


- — / In G 
2-7T1 J r 


t 9 dr 
P( T ) ' 


Additional information follows from the reality of u which leads to u>(X) — u(D) = 
ui(X) — u){D). The reality condition for X implies via Abel’s theorem the existence of a 
meromorphic function R with poles in X + I) and zeros in X + I) (which is a rational 
function in the fundamental polygon). 


(104) 


R(K ) = const 


nti (K - Ej)(K - f~ ) - Qp(KMK) 
ULi( K - Xi){K - Ei) 


where Qo{K) = xo + X\K + ... + xK 9 1 is a polynomial in K with purely imaginary 
coefficients and x = ibe~ 2U . The coefficients x, are related to X via the relation 


(105) (1 


9 9 

-x 2 )f[(K-XMK-Xi) = YliK-EiXK-EJ-QKKKK-OiK-O. 

i =1 i= 1 


We can use the existence of the rational function R to calculate certain integrals of the third 
kind as 


0(u + w(P))0(u + w(P cr ) + cu(£)) 
0(u + u(P a ))Q(u + uj(P) + tu(£)) 


r X+D 


duipo 


R(P) 

P(P CT )' 


(106) 


X+D 
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This makes it possible to give an algebraic expression for S 12 and S 21 . We can write Si ■> 
in the form 
(107) 

_J_ n / [*+ D 

i{£ + £) 00 [Jx+d 


i-D 


dujpQ + / dujpQ 


ID 



duipQ 


with Q independent of P. The second integral can be reexpressed in terms of theta func¬ 
tions and be calculated with the help of so-called root functions: The quotient of two 
theta functions with the same argument but different characteristic is a root function which 
means that its square is a function on Ch ■ Let Pi, i = 1,..., 2g + 2, be the branch points 
of a hyperelliptic Riemann surface Ch of genus g and Aj = w(Pj) with ui(P \) = 0. Fur¬ 
thermore let {*i,..., i g } and {ji ,... ,j g } be two sets of numbers in { 1 , 2 ,..., 2g + 2 }. 
Then the following equality holds for an arbitrary point P £ Ch, 


(108) 


ej^+E Li^K^p)) 

e [Kr + ELi A jk ] ( W (p)) 


= Cl] 


where C\ is a constant independent of K. 
This implies 


(109) exp 

Thus we get with i !106> 

( 110 ) 



Si 2 — 


i(£ + £) \ 2 




l (K-E il )...(K-E iB ) 

{K-E h )...{K-E jg y 


p a 00 + — 


±n 


K-Fj 
K - Ei' 


2—1 


1-x 2 


Ei + Ej) — (£ + £)) + x g ~ 2 


^i=l 


and similarly 

(111) 


<521 


££ 

i{£ + £) 



1 

1 — x 2 




{Ei + Ei) — (£ + £) 



4.3. Counter-rotating dust disk. As an example we will discuss a class of disk solutions 
OTI [32 ] on a genus 2 surface where the disk can be interpreted as two counter-rotating 
components of pressureless matter, so-called dust. The surface energy-momentum tensor 
,S" wi of these models, where a and [3 stand for the t, p and cf> components, is defined on the 
hypersurface £ = 0. The tensor S a P is related to the energy-momentum tensor T n "' which 
appears in the Einstein equations = 8 nTi 11 ' via T“^ = S a/3 e k ~ u 6{(). The tensor 
S' yli can be written in the form 


( 112 ) 


S al3 = cr+u+ul + (J , 


where u± = (1,0, ±fl). We gave an explicit solution for disks with constant angular ve¬ 
locity V. and constant relative density 7 = (<r + — <t_)/(ct++ct_). A physical interpretation 
of the surface energy-momentum tensor S a P will be given in the following section. 

This class of solutions is characterized by two real parameters A and 6 which are related 
to Cl and 7 and the metric potential Uq = U ( 0 , 0 ) at the center of the disk via 


1-7 2 

n 2 


(113) 


A = 2Cl 2 e~ 2U °, S 
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We put the radius po of the disk equal to 1 unless otherwise noted. Since the radius appears 
only in the combinations p/po. (/po and ttpo in the physical quantities, it does not have an 
independent role. It is always possible to use it as a natural length scale unless it tends to 
0 as in the case of the ultrarelativistic limit of the one component disk. The Ernst potential 
will be discussed in dependence of the parameters e = zr /(1 + zr) = 1 — e u ° and 7, 
where zr is the redshift of photons emitted at the center of the disk and detected at infinity. 

The solution is given on a surface of genus 2 where the branch points of the Riemann 
surface are given by the relation E\ = — E 2 and E := E 2 = a + i/3 where 


S fl " ¥ 

(114) a — —1 + —, + 

The function G in < 1941 reads 


(115) 


\J (r 2 - a) 2 + (3 2 + t 2 + 1 

\J (t 2 - a) 2 +/? 2 - (r 2 + 1) 


We note that with a and /3 given, the Riemann surface is completely determined at a given 
point in the spacetime, i.e. for a given value of £. 

Regularity of the solutions in the exterior of the disk restricts the physical parameters to 
0 < 6 < S s (X) := 2 (l + Vl + l/A 2 ) and 0 < A < A c where A c ( 7 ) is the smallest value 
of A for which e = 1. The range of the physical parameters is restricted by the following 
limiting cases: 

Newtonian limit: e = 0 (A = 0), i.e. small velocities Qpo and small redshifts in the disk. 
The function e 2U tends independently of 7 to 1 + A Um, where (/ ; .y is the Maclaurin disk 
solution, and b is of order 9 :! . 

Ultrarelativistic limit : e = 1, i.e. diverging central redshift. For 7 / 1 it is reached 
for A c = 00. The solution describes a disk of finite extension with diverging central 
redshift. For 7=1, the limit is reached for A c = 4.629 .... In this case the solution has a 
singular axis and is not asymptotically flat. This behavior can be interpreted as the limit of 
a vanishing disk radius. With this rescaling the solution in the exterior of the disk can be 
interpreted as the extreme Kerr solution (see EU and references given therein). 

Static limit : 7 = 0 (S = 6 S ( A)). In this limit, the solution belongs to the Morgan and 
Morgan class E) 

One component : 7 = 1 (5 = 0), i.e. no counter-rotating matter in the disk. This is the disk 

of mm. 

Analytic formulas for the complete metric in terms of theta functions are given in 
To evaluate the hyperelliptic integrals in the expressions for the metric we use the numerical 
methods of EU. 

At the disk the branch points £, £ lie on the contour T which implies that care has to 
be taken in the evaluation of the line integrals. The situation is however simplified by the 
equatorial symmetry of the solution which is reflected by the additional involution K —> 
-K of the Riemann surface S 2 for ( = 0. This makes it possible to perform the reduction 
K 2 —> r and to express the metric in terms of elliptic theta functions (see mi)- We denote 
with E w the elliptic Riemann surface defined by p? w = (r 4- p 2 )((r — a) 2 + (3 2 ), and let 
dw be the associated differential of the first kind with u w = i J * In G(-\/r)dw{r). We 
cut the surface in a way that the a-cut is a closed contour in the upper sheet around the cut 
[— p 2 , E\ and that the 6 -cut starts at the cut [00, E\. The Abel map w is defined for P £ T, w 
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as w{P) = dw. Then the real part of the Ernst potential at the disk can be written ; 



(116) 

where 

(117) 


IY 2 ((p 2 + a) 2 + (3 2 ) , , 1 , r _, 

— 2 -2 Y(p 2 + a)+ -^ + dp 


t 

A 7 


Y = 


+ dp 2 fl 2 ju w ) 
y/(p 2 +a) 2 +P 2 $l{u w )' 


It was shown that there exist algebraic relations between the real and imaginary parts of 
the Ernst potential, 


(118) 


<5 2 


(e 4U + b 2 ) = 


— Se 


2U 



- 1 


and the function Z := (a — ao)e 2U : 
(119) Z 2 ~p 2 

Moreover we have 


Se 4U = —e 

A 


2 U 


( 120 ) 


ix 0 = 


Z ( 1/A 


5f \^/ll\ 2 + 5p 2 A 


At the disk, the normal derivatives of the metric functions are discontinuous, but they 
can be expressed in terms of /^-derivatives via 


(e 2C/ )c = 


Z 2 +p 2 +6e AU 


2 Zp 


( 121 ) 


k = ~ 


z 2 


6e 4U 


0 2U 


2 Zp 




This makes it directly possible to determine all quantities in the disk in terms of elliptic 
functions. 

4.4. Metric and Harrison transformation. If stationary axisymmetric solutions for the 
pure vacuum are submitted to a Harrison transformation, the complete transformed metric 
can be constructed. The metric function k is invariant under the action of SU (2,1) trans¬ 
formations. To determine the transformed metric function a', we consider the matrix S in 
©. If we go over from 2 x 2-matrices to 3 x 3-matrices according to the rule 


( 122 ) 


An A12 
A21 A22 


An 0 A12 

0 10 
A21 0 A22 


the matrix Q acts on S as on Thus we get with © 


(123) 


(1 - Q 2 ) 2 


(S12 + 2 iq 2 — q 4 S2i ), 
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which is in accordance with d60l >. This implies with d84l for the function a' 

(124) a'-a' 0 = ~ —-^-(5i2 - g 4 52i). 

(1 - q-y 


To determine a' 0 , one has to consider <Si 2 and S 21 on the axis. In the limit p —> 0, there 
is a non-trivial contribution from the quotient of theta functions in iTTool which diverges as 
1/p. Repeating the considerations of 1181 in the calculation of ao, one finds that the axis 
potentials can be expressed in terms of theta functions on the surface E of genus g — 1 
given by /t 2 = n Li(k — Ei)(K — E.j). Denoting quantities on E with a tilde, one has 
with m on the axis 

(125) ^21 = 0pq(u+2di(oO-)) 

S12 0pq(u + 2cu(oo+)) 


In the equatorially symmetric case, this quotient is identical to one since 2w(oo + ) is a half 
period on E (see ED). Thus we have 


(126) 


a 0 = a 0 


1 + g 4 

(1 - Q 2 ) 2 ' 


To illustrate the class of Harrison transformed hyperelliptic solutions, we will now study 
the transformed counter-rotating dust disks. The metric function f in ( 149b is proportional 
to /, which implies that the transformed solution vanishes exactly where the original so¬ 
lution has zeros. Since the set of zeros of f just defines the ergoregions, the transformed 
solution has the same ergoregions (if any) as the original solution. For the ergoregions of 
the counter-rotating dust disks see El 

For small q, the functions / and b are essentially unchanged since they are quadratic in 
q. The electromagnetic potential $ is in this limit with ED of the form 

(127) $' = —9(1 -f + ib). 


For larger q , /' becomes smaller near the origin. Since its asymptotic values are not 
changed, the growth rate towards infinity increases which is reflected by the mass formula 
El- In the singular limit q —> 1, the function /' is zero for all finite values of z, but 
one at infinity. The behavior for b' is similar with the exception that b' is odd and zero at 
infinity. The function a also becomes singular in the limit q —> 1 which is reflected by the 
diverging factor 1/(1 — q 2 ) 2 and the constant ag d 1 261 which just implies that one can no 
longer choose a to be zero on the axis. In the metric function g' 03 = —a'f, the factors 
(1 — q 2 ) 2 just cancel and the function is only marginally changed with increasing q. The 
typical behavior of g’ 03 for values of q with 0 < \q\ < 1 can be seen in Fig. El The metric 
function is an even function in ( which vanishes on the axis and at infinity. It is analytic 
except at the disk where the normal derivatives have a jump. 

The electromagnetic potential tends to —1 in the limit q —> 1 for finite /, but is zero at 
infinity. The imaginary part is directly proportional to b' as can be seen from m and ED - 
We show a typical situation for values of q with 0 < |g| < 1 in Fig. EJfor the real part and 
in Fig. Elfor the imaginary part. The real part of *1) is an even function in £ which vanishes 
at infinity and has discontinuous normal derivatives at the disk. The imaginary part of $ is 
an odd function in ( and has a jump at the disk. It vanishes at infinity. 

Since /' has the same zeros as /, the transformed solution has a diverging central red- 
shift if the untransformed has, i.e. the ultrarelativistic limits coincide. In the case 7 / 1, 
one has a charged disk of finite extension with diverging central redshift. For 7 = 1 , the 
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FIGURE 4. Real part of the electromagnetic potential $ for e = 0.85, 
7 = 0.95 and q = 0.6. 


solution in the exterior of the disk can be interpreted as an extreme Kerr-Newman metric 
which is obtained as a Harrison-transformed extreme Kerr metric. 
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Figure 5. Imaginary part of the electromagnetic potential <l» for e = 
0.85, 7 = 0.95 and q = 0.6. 


In the Newtonian limit A —> 0, one has / = 1 + XU n and b = A 2 b in lowest order. This 
implies with ( 1491 ) to ( 15 11 1 for 1 — q 2 A 

( 128 ) f = \ + l±*-\u N , b' = ^-^xU, & = - — *—(\U N + i\U), 

1 — q z 1 — q z 1 — q z 

The transformed solution thus has the same Newtonian and post-Newtonian behavior as 
the original metric and in addition an electromagnetic field. The magnetic field is of order 
fl 3 as b'. 

Since the mass is of order A in the Newtonian limit, it is possible to have an extreme 
limit here with M' = Q' as in the Reissner-Nordstrom solution. If we put 1 — q 2 = k\ 
with k > 0, we get in the limit A —> 0 for <1491 and <EB 


(129) 


/' 


(k - Un) 2 ’ 




Un 

k — Un 1 


a static solution similar to the extreme Reissner-Nordstrom solution, but with a jump in the 
normal derivatives of the metric functions at the disk and non-vanishing /' at the origin. 
Since Un < 0 in the whole spacetime, the solution is regular in the exterior of the disk. 
Thus one gets a non-singular limit in the exterior of the disk for q —>• 1 in this case. 

In the static limit one has b' = Ts f l> / = 0, since both are proportional to b. 


(130) 


S' 


(1 -q 2 sr 


<&' = -q 


1 -/ 

1 -q 2 f 


The Harrison-transformed static solution is thus again static with vanishing magnetic field 
but non-zero electric field. 
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5. Energy-momentum tensor 

In this section we will study the energy-momentum tensor at the disk which is a sur¬ 
face layer. For the ‘cut and glue’ techniques, which lead to infinite disks, the energy- 
momentum tensor can be determined using Israel’s invariant junction conditions. In the 
Einstein-Maxwell case we consider the Harrison-transformed counter-rotating disks as an 
example and check the energy conditions. Following GB, the matter is interpreted as 
streams of counter-rotating electro-dust where possible. 

5.1. Disk and energy-momentum tensor. To treat relativistic dust disk, it seems best to 
use Israel’s invariant junction conditions for matching spacetimes across non-null hyper¬ 
surfaces <331 The disk is placed in the equatorial plane and the regions V ± (±£ > 0) are 
matched at the hyperplane £ = 0. This is possible in Weyl coordinates since we are only 
considering dust i.e. vanishing radial stresses in the disk. The spacelike unit normal vector 
of this hypersurface in V + is (n a ) = e k ~ u (0,0,1, 0). The extrinsic curvature Kab of this 
plane in V + is defined as Kab = nA\\B', here capital indices take the values 0, 1, 3 corre¬ 
sponding to the coordinates t, p, £>, | j denotes the covariant derivative with respect to sab, 
the metric on the hypersurface. According to Israel 1331 the jump -/a is = K\ b — K^ B in 
the extrinsic curvature Kab of the hypersurface £ = 0 with respect to its embeddings into 
V ± = {±£ > 0 } is related to the energy momentum tensor Sab of the disk via 

(131) —&'kSab = 1ab — sab1c- 

As a consequence of the field equations the energy momentum tensor is divergence free, 

S AB n ~ 0. 

The relations (11311 lead to 

- Ane^-^Soo = (fc c - 2U C ) e 2U , 

U \Sq 3 — aSoo) = ^-^ a c e2U ■> 

(132) -47re( fc -^(5 3 3-2a5’o3+a 2 Soo) = -~k c p 2 e~ 2U . 

With these formulas it is straight forward to calculate the energy-momentum tensor for a 
given spacetime. The discontinuity of the normal derivatives in the equatorial plane due 
to the ‘cut and glue’ techniques lead as in the Newtonian case to a disk like surface layer. 
There are purely azimuthal tensions in this case. 

Whenever an energy-momentum tensor is worked out in the above way by entering 
with a metric into the Einstein tensor, the question has to be addressed whether the matter 
is physically acceptable. The usual criterion one has to check are the energy conditions, 
see | 66 |. The weak energy condition implies that the energy density must be positive for 
all observers, SabV a V b > 0 where Va is an arbitrary timelike vector. The dominant 
energy condition is satisfied if the weak energy condition holds and if in addition the flux 
of energy is positive for any observer, i.e. that S ab Va is a non-spacelike vector for an 
arbitrary timelike vector Va- An energy-momentum tensor satisfying the dominant energy 
condition will be called physically acceptable, otherwise the matter is exotic (negative 
energy densities, superluminal velocities,...), as Bondi called it, ‘not the cheapest material 
to be bought in the shops’. 

A convenient way to describe the matter of the disk in the pure vacuum case was in¬ 
troduced by Bicak and Ledvinka in 0: An energy-momentum tensor S AB with three 
independent components can always be written as 

S AB = a* p V A V B + p* p W A W B , 


(133) 
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where V and W are the unit timelike respectively spacelike vectors ( V A ) = Ni(l, 0, ui^) 
and where ( W A ) = N 2 (k, 0,1). This corresponds to the introduction of observers (called 
^-isotropic observers (FIOs) in |9| ) for which the energy-momentum tensor is diagonal. 
The condition WaV A = 0 determines k and ui,^ in terms of the metric, 

(134) 

533^00 ~ 500-533 + \/(ff33S00 ~ SOOS 33) 2 + 4(5q3S00 - 500^03) (503-533 ~ 533-S03) 


W0 = 

and 

(135) 


2 (ffo 3 S 33 — 533S03) 


K = — 


5 03 + ^533 


500 + 03 

In non-static cases the FIOs rotate with respect to the locally non-rotating observers for 
whom the metric is diagonal. The latter rotate with the angular velocity u>i with respect to 
infinity 

„„4tl 


(136) 


— _503 _ 
533 


p 2 — a 2 e 4U 


This quantity is a measure for the dragging of the inertial frames with respect to infinity 
due to the rotating matter in the disk. 

If the dominant energy-condition holds and if the pressure is positive, one has Pp/cr* < 
1, and the matter in the disk can be interpreted as in ED either as having a purely azimuthal 
pressure or as being made up of two counter-rotating streams of pressureless matter with 
proper surface energy density a* /2 which are counter-rotating with the same angular ve¬ 
locity y/ Pp/cr* (which is below 1, the velocity of the light), 


(137) 


S AB = ^a*(U A U B + U A U B ) 


where (U±) = U*(v A ± s Jp*/a*w A ) is a unit timelike vector. We will always adopt the 

latter interpretation if the condition pl/c t* < 1 is satisfied which is the case in the example 
ED The energy-momentum tensor <11371 is just the sum of two energy-momentum tensors 
for dust. Furthermore it can be shown that the vectors U± are geodesic vectors with respect 
to the inner geometry of the disk: this is a consequence of the equation S ab \\b = 0 
together with the fact that U± is a linear combination of the Killing vectors. 

Using the ‘cut and glue’ formalism, disk sources for all known stationary axisymmet- 
ric metrics can be given. Sources for static spacetimes were discussed in (7), the Kerr 
spacetime was considered in Ell- If the strip cut-off the spacetime includes the horizon 
completely in the latter case, the matter in the disk can be interpreted as dust. 

In the presence of electromagnetic fields, a discontinuous electromagnetic tensor F a/3 
leads to a current density J a via the Maxwell equations 

(138) F A f *. 4 3 = ~^= {V=gF A P) >f) = -47 tJ a . 

Contributions to J A arise only via the normal derivatives at the disk. We define the current 
density j A in the disk as s AB by the relation J A =: j A S(Q. The electromagnetic 

energy-momentum tensor does not produce a ci-type contribution to the Einstein-Maxwell 
equations since F a P is bounded at the disk. With m we get using equatorial symmetry 
(the derivatives are taken at £ = 0 + ) 

27rjo = -(5R$)c, 2 tt (j 3 - a jo) = -y(3$) p . 


(139) 
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The continuity equation (the Bianchi identity) at the disk = F 1 "'= 0 leads to 

the condition 

(140) 5 oo, P s 00 + 2503 , P s 03 + g 33 , p s 33 = 2 (F 10 j° + F 13 j 3 ). 

We remark that one can substitute one of the equations (II 321) by (11401 in the same way 
as one replaces one of the field equations by the covariant conservation of the energy 
momentum tensor in the case of three-dimensional perfect fluids. This makes it possible to 
eliminate from m and to treat the energy-momentum tensor at the disk purely on the 
level of the Ernst equation. It is straight forward to check the consistency of this approach 
with the help of (ED- Thus one can solve boundary value problems of dust disks without 
using k. 

The interpretation of the energy-momentum tensor is not changed with respect to the 
pure vacuum case. However, free particles will no longer move on geodesics but on so- 
called electro-geodesics where the Lorentz force is added to the geodesic equation which 
leads to for each component of dust. It follows that the matter streams will only 
move on electro-geodesics in the FIO frame if jAW A = 0, i.e. if there are no currents in 
this frame. This is in general not the case which implies that the FIOs cannot interpret 
the matter in the disk as freely moving charged particles even if the energy conditions are 
satisfied which was always possible in the pure vacuum case. However the matter can still 
be interpreted as a fluid with a purely azimuthal pressure. 

Since the splitting of the matter into two streams of dust is not unique, an alternative 
approach is to interpret the matter in the disk as two streams of particles moving on electro¬ 
geodesics in the asymptotically non-rotating frame. To this end we make the ansatz Co) 


(141) s AB = a+U A U B + a~U A U B , j A = atU A + a~U A 

with (U±) = N±( 1,0, uj±). The angular velocity follows from the electro-geodesic equa¬ 
tion for each component, 

(142) ^ a m(9o 0,p + 2503, p w± + 533, pU±) = cr^(H 0 ,p + A 3 ^pU)±). 

Such an interpretation is possible if the angular velocities are real. Moreover the velocities 
ui±p in the disk should be smaller than 1 to avoid superluminal velocities, and the energy 
densities have to be positive, which means we are not interested in tachyonic or otherwise 
exotic matter. Relation QD leads to 


(143) 


crtNl = 


UJ-\- — UJ- 


N 2 = 

UJ- — CJ-\- 


t+N 2 + = 


CJ-\- — UJ- 


(144) 
and 

(145) 

If we enter (fT42t with this, we obtain 


s 00 

a ~N 2 = - _ + 

u m iy — ‘ 

UJ- — UJ-\- 


s 33 - cu+s 03 
s 03 — cu+s 00 


-T 2 ± JTI - 7’| 7 3 
Ti 


(146) 


w± = 
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with 


,0 „03 _ „-3 „00 

O A J * J ° 

y33,p Z/1 3,p s 03 s 03 _ s oo s 33 

j°S 03 — j 3 S°° j 3 S 03 - j°S 33 

903,p — j4 °,P g 03 s 03 _ s 00 s 33 — ^ 3 ’ P s 03 s 03 — S 00 S 33 
~3 „03 _ -0 q 33 

g00,p zyi 0 ,P s 03 s 03 _ s 00 s 33 ‘ 

The densities then follow from ED where the continuity equation (TT40I guarantees that 
this system can be solved. 

Using the above formalism, disk sources for conformastationary spacetimes ED, mag¬ 
netostatic metrics QU and the Kerr-Newman family El have been constructed and dis¬ 
cussed. 


T\ = 
T 2 = 
(147) T 3 = 


5.2. Energy-momentum tensor of the Harrison-transformed counter-rotating dust 
disk. As an example we will now discuss the Harrison-transformed counter-rotating dust 
disk. It remains to be checked whether and in which range of the parameters the energy- 
momentum tensor fl50l is physically acceptable. The above discussion of the metric in¬ 
dicates the extreme behavior of the metric functions for q close to one. It is plausible that 
the matter in the disk which is in the present example the source of such an extreme metric 
will in general not be physically acceptable. There can be maximal q smaller than 1 for 
given A and 5 which limits the physical range of the parameters. 

To discuss the energy-momentum tensor and the currents in the disk, it is helpful to 
use the algebraic relations ( II 1 8b to (fl20t between the metric functions which exist at the 
disk, and which imply similar relations between the transformed potentials. With (II 1 01 
and ED, we get for S21 


(148) 


S21 = ££S 12 +ix 0 f, 


and thus with ED for the metric function a' 
(149) 

(1 -q 2 ) 2 (a-a 0 y = {a - a 0 ) l + g 4 | | 


1/A 2 + S 
Vl/A 2 + 5p 2 




where a' Q is given by ED- With this function we can can calculate the angular velocity 
of the locally non-rotating observers u>i ED - The dependence of uji on e and 7 has been 
discussed in ED- As a function of q it is monotonically decreasing as can be seen in 
Fig- ID The reason for this behavior is that the function /' tends to zero in the limit q —> 1 
for finite p, £ whereas g' 03 changes shape but remains finite. Thus the overall behavior 
of u>i is dominated by /'. The deformation of the function (j[ yi via q has, however, the 
consequence that 07 has its maximum for large q no longer at the center at the disk but near 
the rim. 

The energy-momentum tensor at the Harrison-transformed disk can be calculated via 
d 1 32b . Expressing the right-hand sides with the help of d49b and <1501 via the original 
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Figure 6 . Angular velocity of the locally non-rotating observers with 
respect to infinity for e = 0.85, 7 = 0.95, and <7 = 0 ,..., 0.9. 


functions, we get 


(150) 4 

where N = (1 - 
This implies 


= { - L ^(^b p -(l- q i f+ q *b 2 )f c +2 q *fbb^ 

-- ((! - 9 4 / 2 + Q*b 2 )b p + 2q A bff p ) , 

2(1 -q*ypz p ' 

9 2 /) 2 + q 4 b 2 . 


Znjo = q ^ 1 N 2 ) (-((! - g 2 /) 2 - q 4 b 2 )f( + 2 g 2 h(l - q 2 .f)b c ) 

(151^7t(j 3 - ajo) = ij—^j^j(2q 2 b(l-q 2 f)f p + ((l~q 2 f) 2 -q 4 b 2 )b p ). 

We will interpret the matter in the disk as in Q 2 l in two ways: The FIOs for whom the 
tensor s AB is diagonal rotate with angular velocity with respect to infinity. We show 
uj 9 for several values of q in Fig.0 It can be seen that oj 9 decreases monotonically with 
q. For large q the maximum of the angular velocity is near or at the rim of the disk in this 
example. 

Since the energy conditions are satisfied here, the FIOs can interpret the matter in the 
disk as a fluid with a purely azimuthal pressure. Alternatively they can interpret it as being 
made up of two streams of counter-rotating pressureless matter, where the velocity f \ c p in 
the streams is below the speed of light. However it can be seen in Fig. [ 8 ] that there are in 
general currents in the frame of the FIOs. Thus an interpretation of the matter by the FIOs 
as freely moving charged particles is not possible. 
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Figure 7. Angular velocity of the FIOs with respect to infinity for e = 
0.85 ,7 = 0.95, and q = 0, .., 0.9. 



Figure 8 . Currents in the frame of the FIOs for e = 0.85, 7 = 0.95, 
and <7 = 0 ,..., 0.9. 


The second interpretation in terms of two streams of electro-dust is limited by the con¬ 
ditions that the angular velocities have to be real and that the energy densities have to be 
positive. Numerically one finds in the present example that the angular velocities are real. 
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but in a wide range of the parameters there are negative energy densities and tachyonic 
behavior. Already in the uncharged case there are infinite velocities in strongly relativistic 
settings with negligible counter-rotation which are due to extrema of the metric function 
<733 in the disk. In this case the quantity T\ in (fim is zero which leads to a diverging cu_. 
Increasing q only enhances this effect. The result is that an interpretation as non-tachyonic 
counter-rotating matter on electro-geodesics with positive energy densities is only possible 
if q, e and 7 are not too large. In other words large values of q are in this setting only 
possible in post-Newtonian or nearly static situations. We show plots of the angular ve¬ 
locities uj± in Fig.[9]for e = 0.36 and 7 = .08 where values of q up to 0.75 are possible. 
The corresponding densities are given in Fig. HSI and the charge densities in Fig. 1111 The 


s" 0 - 


q=0 


- q=0.g ~ 

q=0.4 


q-0.6 


q=0.75 


q=0.75 

q=0.6 


q-0.4 


-q=0.2 

q=0 


0.1 0.2 0.3 0.4 0.5 0.6 

P 


Figure 9 . Angular velocities oj± for e = 0.36, 7 = 0.08, and <7 = 0 ,0.2,..., 0.8. 

densities vanish always at the rim of the disk, the charge densities are identically zero for 
<7 = 0. For larger values of q, the angular velocity becomes bigger and bigger at the 
rim of the disk and finally diverges. The density < 77 , is in this case negative in the vicinity 
of the rim of the disk. 

6. Einstein-Maxwell equations and the Riemann-Hilbert problem 

In this section we present recent results on theta functional solutions to the Einstein- 
Maxwell equations which are constructed with Riemann-Hilbert techniques. We first ob¬ 
tain the hyperelliptic solutions starting from so-called Schlesinger equations. This method 
is then extended to the Einstein-Maxwell case. The keypoint in the construction is the 
solution of the Riemann-Hilbert problem on multi-sheeted Riemann surfaces in terms of 
Szego kernels by Korotkin G3- The material in this section is based on E2. 

6.1. Riemann-Hilbert problem and Schlesinger equations. Theorem 3.1 in section 3 
offers the possibility to construct matrices with certain analytic properties in the spectral 
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Figure 10. Matter densities <7+ and cr m (dashed) for e = 0.36, 7 = 
0.08, and q = 0,0.4,0.8. 



Figure 11. Charge densities cr+ and a e (dashed) for e = 0.36, 7 = 
0.08, and q = 0, 0.4,0.8. 


parameter which lead to solutions to the Einstein-Maxwell equations. The task is to solve a 
so-called Riemann-Hilbert problem for T to obtain a matrix with prescribed singularities. 

The Riemann-Hilbert problem, also known as Hilbert’s 21st problem, can be stated as 
follows: Consider a linear Fuchsian ordinary differential system on CPI for some M x M 
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matrix T, 


(152) 


d'y 


Ai 


N 

V- 

i=i 7 — Ij 




where the matrices Aj £ sl(M, C) are independent of 7 , with initial condition ^(7 = 
00 ) = I. At the singularities of the equation, the solutions will not be single-valued on 
CPI. If one goes on a small loop around the singularities, a matrix '1 which is a fundamen¬ 
tal solution of the differential system will change by multiplication with some matrix which 
is called the monodromy matrix. The Riemann-Hilbert problem (inverse monodromy prob¬ 
lem) is the construction of the differential equation for a given monodromy matrix, for 
more details see 1681 . 

The mathematical formulation of this problem for an M x M matrix reads: fix an 
SL(M, C) monodromy representation Ad of the fundamental group tti[CP 1 \{ 7 i, ... , 7 at}]- 
Choose the standard set of generators Zi,..., Zjv of 7 Ti [CPI \ { 71 ,..., 7 at}] such that the 
contour l rn encircles only one singularity 7 m , and the relation / M ... l x - I is satisfied. 
Then the monodromy representation M is defined by the corresponding set of N SL(M) 
matrices Adi,..., Adjy satisfying the relation Adjv • ■ • Adi = I. The Riemann-Hilbert 
problem is to find a SL(M, C)-valued function ’I'(P) defined on the universal covering X 
of CPI \ { 71 ,..., 7 jv} such that'T gains the right multiplier equal to Ad m being analyti¬ 
cally continued along contour l m . 

Note that in this formulation the solution of the Riemann-Hilbert problem is not unique; 
different solutions are related by so-called Schlesinger transformations (multiplications 
from the left with appropriate matrix-valued rational functions). It is known that the 
Riemann-Hilbert problem cannot be solved explicitly in terms of known special functions 
in the generic case. The largest class of explicitly solvable problems known so far is in the 
class of problems with quasi-permutation monodromies (when each row and each column 
of monodromy matrix contain exactly one non-vanishing entry). The universal cover X is 
in this case a M-sheeted compact Riemann surface. 

A solution T £ SUM. C) solves the equation (11521 . The behavior near the singulari¬ 
ties is given (in general position, i.e. the difference between the eigenvalues is non-integer 
for each of the matrices Aj) by the asymptotic expansion of > 1 , 

(153) ^( 7 ) = Qj(I + 0(7 - - Ijf’Cj, 

where Qj,Cj £ SL(M,C) and where Tj is a diagonal, tracefree matrix. The matrices 
Mj = Cj e 27 r, ' T > C 3 are the monodromy matrices. The function T is thus only single¬ 
valued on the universal covering of CPl\{ 7 i,..., jn }• The so-called isomonodromy 
condition is satisfied if the monodromy matrices are independent of the ■y ] which implies 
that one can change the position of the singularities in (11521 without affecting the mon¬ 
odromy behavior. In general position, this condition implies 

(154) ^ = 

«7i 7-7) 

The consistency of these equations with <n~52l leads to the Schlesinger system |[69;| for the 
residues Ai, 


(155) d X = 


[Ai , Aj ] 7/7 dAi 

) * r Ji —] — — 

7 i - Ij aji 


j¥* 


[Aj, Aj] 

7 i - Ij ' 


To solve the Riemann-Hilbert problem for quasi-permutation monodromies 1191 . we in¬ 
troduce two further objects on X: the prime form P(P, Q) has the property that it vanishes 
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at exactly one point on X, namely when P = Q. Since the surface is compact, there can be 
no function with this property (the Riemann vanishing theorem states that Riemann theta 
functions have exactly g zeros). The prime form is a (— —^) differential on the bundle 

X x X 


(156) 


E(a, b) 


o*(L b ) 

h&(a)h/±(b) 


where * = [p*q*] denotes a non-singular odd characteristic, and where the spinor is 
given by h 2 ^{a) = Y2a=i dz a ®*(®)d ul a( T a) (r Q denotes local coordinates in the vicinity 
of a point a). In local coordinates the prime form reads 


E(a, b) 


T a - n 

Vdi^Vdn 


The Szego kernel ,S pq (7 ;) . Q) is defined as 


(157) 


5 pq (P,Q) 


9pq( z + fp ) 

0 pq (z )E(P,Q) 



rp - tq 


As can be seen from the representation in local coordinates, it can be viewed as a general¬ 
ization of the Cauchy kernel to Riemann surfaces. 

In terms of the Szego kernel the solution to the Riemann-Hilbert problem was written 
by Korotkin ED in the form 

(158) $ 3i (H,A)=S Pq (A^,A W )£o(H,A) = 


where K 1 -- 11 denotes a point on the jth sheet with projection K on CPI, and where 

K-X 


(159) 


Eo(K,X) = 


\fdK\fd\ 

is the prime form on CPI. Obviously 7 defined in im is not a single-valued function on 
CPI (it is only single-valued on X') since it is a solution of the Riemann-Hilbert problem. 


6.2. Vacuum case. It is instructive to reconsider first the two-dimensional vacuum case. 
In section 4 we treated the Ernst equation as the integrability condition for the linear system 
d79> which had the advantage that the Ernst potential is given as one component of T(oc + ). 
Several linear systems for the Ernst equation are known in the literature which are related 
through gauge transformations, see 1701 . In this section we will use a linear system which 
makes the symmetry of the Ernst equation obvious since the matrix of the linear system is 
an element of the coset space SL(2,M)/SO(2) (in an abuse of notation we will call this 
matrix also T), 

g.g - 1 g^g - 1 

(160) , % = -f— 4s 

1 — 7 5 1 + 7 

where 


(161) 7 (K, 0 = + , 

and where 


(162) 


g 


1 

£ + £ 


2 _ i{£-_£) \ 

i{£ - £) 2££ ) ' 


It is straight forward to adapt theorem 3.1 to the linear system J 1 601 and to the solution of 
a Riemann-Hilbert problem as is done in the following theorem by Korotkin and Nicolai 
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ED: 


Theorem 6.1. Let the matrices Aj £ sl( 2,C) satisfy the Schlesinger system (1/551 and 
let 4/ ( 7 ) be the corresponding solution of \152\ . Suppose that the matrix 4/ satisfies the 
additional conditions 

(163) «' t (1/7)«'^ 1 (0)^(7) = I 
and 

(164) (- 7 ) = 4 /( 7 ). 

Let "/j = 7 (Kj,£,£) with Kj £ C independent off £. Then the matrix 

(165) Q = $(7 = 0,^,0 

is real and symmetric and satisfies the Ernst equation, and the function satisfies the 
system ( 17601 ). 

The proof is analogous to the one for theorem 3.1, see EH The solution of the 
Riemann-Hilbert problem gives a 57. (2. C) matrix 47, the conditions d 1 631 and iTT64l en¬ 
sure that It is in the coset SL(2,M.)/SO(2). The involution a which interchanges the 
sheets on C acts as 7 —> 1 / 7 , condition dl63b is thus equivalent to the reduction condition 

III. The reality condition dl~64b was previously incorporated in the normalization condition 

IV. 

In the two-dimensional vacuum case, the covering surface X is hyperelliptic given by 
p 2 = IX/Li (7 — 7j) (we assume that N is even). Since 7 depends on the spectral parameter 
K via d 161 b which lives on the two-sheeted surface C, the surface is the four-sheeted cov¬ 
ering of the complex plane C with Hurwitz diagram Fig. □of section 4. The factorization 
£/a leads as before to the hyperelliptic surface Ch- 

The corresponding solution of the Ernst equation can be obtained via (11651 One obtains 


(166) 


£ = Qpq(^(oO + )) 

©pqMoO-)) ' 


The relation to the previous form of the Ernst potential can be established as follows: Let 
g = g + n and Pg+j = hj £ R, Qg+j = 0 for j = 1,... ,n. Consider the limit of collapsing 
branch cuts for j > g, i.e. Eg + j —» In this limit all quantities entering d 1 66b can be 

expressed in terms of quantities (denoted with a tilde) on the surface Ch of genus g given 
by p 2 = (K — £)(A' — £) Y\i=i(K ~ Ei){K — Ft), the surfaces Ch with the collapsing 
cuts removed. The holomorphic differentials have the limit, see 1621, 


(167) duti —> duti, i = l,...,g, duti —> dui E - E +, i=g + l,...,n. 

In other words the holomorphic differentials of Ch become holomorphic differentials on 
Ch and differentials of the third kind with poles at the collapsed branch cuts. Since the 
6 -periods of differentials of the third kind can be expressed in terms of the Abel map of the 
poles, see e.g. (72). 

I dujpQ = u>j(P) - u>j(Q), 

b 3 


(168) 
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one can use formula (19Tb to get for (11661 
(169) 

£ _ BpqfuH 00 ^ + E"=i hj{u(E-)~ u>{E+))) ^ (sr^ h 

©pq(^(o°“) + E"=1 h j(u(E~)-w(Et))) P ypl 3 

By taking the limit Y^j=i / r from a sum to a line integral over the Ej, we get after a 
partial integration (we assume In G vanishes as the limits of integration) and the identifica¬ 
tion h{K) = <9ff(ln G) formula (19 3 ) for ( 1 1 69> where we have used f^dujpQ = Jp cIujab- 
For details of the above construction see tZ3. 

6.3. The Einstein-Maxwell case. The construction of algebro-geometric solutions to the 
Einstein-Maxwell equations can be carried out as above for the vacuum case. The associ¬ 
ated linear system has the form d 1 521 where this time the matrix Q £ 5(7(2, l)/5(7(2). 
The matrix 'b is constructed as the solution of a Riemann-Hilbert problem as before. This 
time the covering surface is three-sheeted as noted in section 3. Since the spectral pa¬ 
rameter 7 varies on the two-sheeted surface C , the covering surface £ is now six-sheeted. 
Factorizing with respect to the involution a of C leads to a three-sheeted surface £3 which 
replaces the hyperelliptic surface Ch in the vacuum case. All branch points of this surface 
are constant except the two branch points £ and £ (the fixed points of the involution <r), 
where the first two sheets are glued. If the third sheet ‘detaches’ from the first two, one 
gets the vacuum solutions. 

The reduction conditions to obtain Q £ 5(7(2,1)/5(7(2) imply that £3 is invariant 
with respect to the holomorphic involution < 7 , acting on every sheet of £3 as 7 —> 1 /"f. In 
addition £3 has to be invariant with respect to the anti-holomorphic involution r, acting as 
7 —> —7 on the third sheet of £3; on the first and second sheets r acts as a superposition 
of the conjugation 7 —> —7 with interchange of the first and second sheets. 

The solution of the Riemann-Hilbert problem is again provided by m for M = 3. 
The Ernst potentials can be obtained from ’3/( 00 ^) as before. Details of the construction 
will be published in a forthcoming paper by Korotkin. The class of algebro-geometrical 
solutions of the Einstein-Maxwell system constructed in El is a partial case of the above 
construction when the third sheet is attached to the first and second sheets only via branch 
points of multiplicity two (i.e. all three sheets are glued together at these points). 

7. Conclusion and Outlook 

In this article we reviewed exact solutions to the Einstein-Maxwell equations with disk 
sources. The ‘cut and glue’ techniques start with known exact solution of the equations 
with equatorial symmetry. The hyperplanes Co and —Co are identified. This corresponds 
to removing a strip from the spacetime. At the newly formed equatorial plane, the normal 
derivatives of the metric functions will be discontinuous which leads via the field equa¬ 
tions to a (5-type energy-momentum tensor. These infinitesimally thin disks have infinite 
extension. For asymptotically flat spacetimes, the mass of the disks is however finite. 
The matter in the disk has to satisfy the energy conditions in order to be physically ac¬ 
ceptable. If this is the case, an interpretation of the matter in the disk as two streams of 
counter-rotating dust is possible if the velocities in the disk are subluminal. In the absence 
of electromagnetic fields, the particles of the streams move on geodesics. In the non-static 
Einstein-Maxwell case, the electro-geodesic equations leads to additional conditions on the 
matter which can only be satisfied in a certain range of the physical parameters. ‘Cut and 
glue’-techniques thus provide a generally applicable method to find (infinite) disk sources 
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to known spacetimes. They do not provide, however, a method to construct new solutions 
to the Einstein-Maxwell equations. 

The complete integrability of the stationary axisymmetric Einstein-Maxwell equations 
offers powerful techniques to generate new solutions. The most efficient methods known 
up to day arise from algebraic geometry and lead to explicit solutions in terms of theta 
functions associated to certain Riemann surfaces. In the pure vacuum case, these surfaces 
are hyperelliptic, which makes the use of the powerful hyperelliptic calculus possible. 
These solutions contain the ‘solitonic’ solutions as the Kerr solution as limiting cases. 
An important feature of the algebro-geometric solutions is that general regularity theorems 
can be established. Thus the question of global regularity in the exterior of the disk source 
can be addressed. With these techniques, disks of finite extension can be constructed as 
the presented family of counter-rotating dust disks which leads to new solutions to the 
Einstein equations. The task is then to solve boundary value problems at the disk arising 
from physical models. 

In the Einstein-Maxwell case, hyperelliptic disk solutions were constructed by exploit¬ 
ing the SU (2, l)-symmetry of the field equations. This leads to charged disks. By studying 
the asymptotics of the so generated spacetime, it was shown that the disks always have a 
non-vanishing total charge and a gyromagnetic ratio of 2. To generate solutions without 
total charge but with non-trivial magnetic field, it seems necessary to study theta functions 
on three-sheeted surfaces. We gave a review on recent results on the Riemann-Hilbert 
problem on multi-sheeted surfaces which can be applied to the Einstein-Maxwell case. 

In order to obtain disk solutions on three-sheeted surfaces, the solutions to Riemann- 
Hilbert problems have to be adapted to the Einstein-Maxwell equations. The general struc¬ 
ture of the physically interesting Riemann surfaces has to be explored. An important point 
is to establish the analyticity properties of the solutions. If this is achieved, disk solutions 
can be studied. In order to construct solutions to prescribed matter models, one will have 
to adopt the algebraic approach 13211741 to the three-sheeted case. For a complete under¬ 
standing of the solution, it will be necessary to extend the numerical code for hyperelliptic 
theta functions 1541 to these surfaces. 

I thank D. Korotkin for many discussions and hints, and M. King for carefully reading 
the manuscript. This work was supported by the Schloessmann foundation. 
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